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Abstract. The discovery of X-rays and radioactivity in the waning years of the 19th 
century lead to one of the most awe inspiring scientific eras in human history. The 
20th century witnessed a level of scientific discovery never before seen or imagined. 
At the dawn of the 20th century only two forces of Nature were known - gravity 
and electromagnetism. The atom was believed by chemists to be the elemental, 
indestructible unit of matter, coming in many unexplainably different forms. Yet J.J. 
Thomson, soon after the discovery of X-rays, had measured the charge to mass ratio 
of the electron, demonstrating that this carrier of electric current was ubiquitous and 
fundamental. All electrons could be identified by their unique charge to mass ratio. 

In the 20th century the mystery of the atom was unravelled, the atomic nucleus 
was smashed, and two new forces of Nature were revealed - the weak force [responsible 
for radioactive [3 decay and the nuclear fusion reaction powering the stars] and the 
nuclear force binding the nucleus. Quantum mechanics enabled the understanding of 
the inner structure of the atom, its nucleus and further inward to quarks and gluons 
[the building blocks of the nucleus] and thence outward to an understanding of large 
biological molecules and the unity of chemistry and microbiology. 

Finally the myriad of new fundamental particles, including electrons, quarks, 
photons, neutrinos, etc. and the three fundamental forces - electromagnetism, the 
weak and the strong nuclear force - found a unity of description in terms of relativistic 
quantum field theory. These three forces of Nature can be shown to be a consequence 
of symmetry rotations in internal spaces and the particular interactions of each particle 
are solely determined by their symmetry charge. This unifying structure, describing 
all the present experimental observations, is known as the standard model. Moreover, 
Einstein's theory of gravity can be shown to be a consequence of the symmetry of local 
translations and Lorentz transformations. 

As early as the 1970s, it became apparent that two new symmetries, a grand 
unified theory of the strong, weak and electromagnetic interactions in conjunction 
with supersymmetry, might unify all the known forces and particles into one unique 
structure. Now 30 years later, at the dawn of a new century, experiments are on 
the verge of discovering (or ruling out) these possible new symmetries of Nature. 
In this article we try to clarify why supersymmetry [fi] and supersymmetric grand 
unified theories [SUSY GUTs ] are the new standard model of particle physics, i.e. the 
standard by which all other theories and experiments are measured. 
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1. Introduction 

Supersymmetry is a space-time symmetry; an extension of the group of transformations 
known as the Poincare group including space-time translations, spatial rotations 
and pure Lorentz transformations. The Poincare transformations act on the three 
space x and one time coordinate t. The supersymmetric extension adds two anti- 
commuting complex coordinates 9 a , a — 1,2 satisfying 9 a 9p + 9p9 a = 0. Together 
they make superspace z = {t,x, 9 a } and supersymmtry transformations describe 
translations/rotations in superspace. Local supersymmetry implies a supersymmetrized 
version of Einstein's gravity known as supergravity. 

In the standard model, ordinary matter is made of quarks and electrons. All of 
these particles are Fermions with spin s = satisfying the Pauli exclusion principle. 
Hence no two identical matter particles can occupy the same space at the same time. 
In field theory, they are represented by anti-commuting space-time fields, we generically 
denote by, ip a (x,t). On the other hand, all the force particles, such as photons, 
gluons, W ± , Z°, are so-called gauge Bosons with spin s = lh, satisfying Bose-Einstein 
statistics and represented by commuting fields <f>(x,t). As a result Bosons prefer to sit, 
one on top of the other; thus enabling them to form macroscopic classical fields. A 
Boson - Fermion pair form a supermultiplet which can be represented by a superfield 
*&(z) = <f>(x, t)+9 ip(x, t). Hence a rotation in superspace, rotates Bosons (force particles) 
into Fermions (matter particles) and vice versa. 

This simple extension of ordinary space into two infinitesimal directions has almost 
miraculous consequences making it one of the most studied possible extensions of the 
standard model [SM] of particle physics. It provides a "technical" solution to the so- 
called gauge hierarchy problem, i.e. why is Mz/M p i <C 1. In the SM, all matter derives 
its mass from the vacuum expectation value [VEV] v of the Higgs field. The W ± and 
Z° mass are of order g v where g is the coupling constant of the weak force. While 
quarks and leptons (the collective name for electrons, electron neutrinos and similar 
particles having no strong interactions) obtain mass of order A v where A is called a 
Yukawa coupling; a measure of the strength of the interaction between the Fermion and 
Higgs fields. The Higgs vacuum expectation value is fixed by the Higgs potential and in 
particular by its mass uih with v ~ mu- The problem is that in quantum field theory, 
the Lagrangian (or bare) mass of a particle is subject to quantum corrections. Moreover 
for Bosons, these corrections are typically large. This was already pointed out in the 
formative years of quantum field theory by [Weisskopf (1939)]. In particular for the Higgs 
we have m 2 H = m§+a A 2 where itlq is the bare mass of the Higgs, a represents some small 
coupling constant and A is typically the largest mass in the theory. In electrodynamics a 
is the fine-structure constant and A is the physical cutoff scale, i.e. the mass scale where 
new particles and their new interactions become relevant. For example, it is known that 
gravitational interactions become strong at the Planck scale M pt ~ 10 19 GeV; hence we 
take A ~ M p i. In order to have M z <C M p i the bare mass must be fine-tuned to one 
part in 10 17 , order by order in perturbation theory, against the radiative corrections in 
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order to preserve this hierarchy. This appears to be a particularly "unnatural" accident 
or, as most theorists believe, an indication that the SM is incomplete. Note that neither 
Fermions nor gauge Bosons have this problem. This is because their mass corrections 
are controlled by symmetries. For Fermions these chiral symmetries become exact 
only when the Fermion mass vanishes. Moreover with an exact chiral symmetry the 
radiative corrections to the Fermion's mass vanish to all orders in perturbation theory. 
As a consequence when chiral symmetry is broken the Fermion mass corrections are 
necessarily proportional to the bare mass. Hence mp = m +a m log(A/mo) and a light 
Fermion mass does not require any "unnatural" fine-tuning. Similarly for gauge Bosons, 
the local gauge symmetry prevents any non-zero corrections to the gauge boson mass. 
As a consequence, massless gauge bosons remain massless to all orders in perturbation 
theory. What can we expect in a supersymmetric theory? Since supersymmetry unifies 
Bosons and Fermions, the radiative mass corrections of the Bosons are controlled by the 
chiral symmetries of their Fermionic superpartners. Moreover for every known Fermion 
with spin \h we necessarily have a spin Boson (or Lorentz scalar) and for every spin 
lh gauge Boson, we have a spin gauge Fermion (or gaugino). Exact supersymmetry 
then requires Boson- Fermion superpartners to have identical mass. Thus in Ban electron 
necessarily has a spin superpartner, a scalar electron, with the same mass. Is this a 
problem? The answer is yes, since the interaction of the scalar electron with all SM 
particles is determined by B. In fact, the scalar electron necessarily has the same charge 
as the electron under all SM local gauge symmetries. Thus it has the same electric 
charge and it would have been observed long ago. We thus realize that Bean only be 
an approximate symmetry of Nature. Moreover it must be broken in such a way to 
raise the mass of the scalar partners of all SM Fermions and the gaugino partners of 
all the gauge Bosons. This may seem like a tall order. But what would we expect to 
occur once Bis softly broken at a scale A fi ? Then scalars are no longer protected by the 
chiral symmetries of their Fermionic partners. As a consequence they receive radiative 
corrections to their mass of order 5m 2 oc a Ag log(A/m ). As long as A B < 100 TeV, 
the Higgs Boson can remain naturally light. In addition, the gauge Boson masses are 
still protected by gauge symmetries. The gauginos are special, however, since even if 
Bis broken, gaugino masses may still be protected by a chiral symmetry known as R 
symmetry [Farrar and Fayet (1979)]. Thus gaugino masses are controlled by both the 
Band R symmetry breaking scales. 

Before we discuss SUSY theories further, let us first review the standard model 
[SM]in some more detail. The standard model of particle physics is defined almost 
completely in terms of its symmetry and the charges (or transformation properties) of 
the particles under this symmetry. In particular the symmetry of the standard model 
is SU(3) x SU(2) x U(1)y- It is a local, internal symmetry, by which we mean it 
acts on internal properties of states as a rotation by an amount which depends on the 
particular space-time point. Local symmetries demand the existence of gauge Bosons 
(or spin 1 force particles) such as the gluons of the strong SU (3) interactions or the 
W ± , Z° or photon (7) of the electroweak interactions SU(2) x U(l) Y - The strength of 
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the interactions are determined by parameters called coupling constants. The values of 
these coupling constants however are not determined by the theory, but must be fixed 
by experiment. 

There are three families of matter particles, spin 1/2 quarks and leptons; each 
family carrying identical SM symmetry charges. The first and lightest family contains 
the up (u) and down (d) quarks, the electron (e) and the electron neutrino (v e ) (the 
latter two are leptons). Two up quarks and one down quark bind via gluon exchange 
forces to make a proton, while one up and two down quarks make a neutron. Together 
different numbers of protons, neutrons bind via residual gluon and quark exchange forces 
to make nuclei and finally nuclei and electrons bind via electromagnetic forces (photon 
exchanges) to make atoms, molecules and us. The strong forces are responsible for 
nuclear interactions. The weak forces on the other hand are responsible for nuclear 
beta decay. In this process typically a neutron decays thereby changing into a proton, 
electron and electron neutrino. This is so-called f3~ decay since the electron (or beta 
particle) has negative charge, — e. [3 + decays also occur where a proton (bound in the 
nucleus of an atom) decays into a neutron, anti-electron and electron neutrino. The 
anti-electron (or positron) has positive charge, +e but identical mass to the electron. If 
the particle and anti-particle meet they annihilate, or disappear completely, converting 
their mass into pure energy in the form of two photons. The energy of the two photons 
is equal to the energy of the particle - anti-particle pair, which includes the rest mass 
of both. Nuclear fusion reactions where two protons combine to form deuterium (a p-n 
bound state), e + + v e is the energy source for stars like our sun and the energy source 
of the future on earth. The weak forces occur very rarely because they require the 
exchange of the W ± , Z° which are one hundred times more massive than the proton or 
neutron. 

The members of the third family {t, b, r, u r } are heavier than the second family 
{c, s, /i, z/ M } which are heavier than the first family members {u, d, e, v e }. Why 
there are three copies of families and why they have the apparent hierarchy of masses 
is a mystery of the SM. In addition why each family has the following observed charges 
is also a mystery. A brief word about the notation. Quarks and leptons have four 
degrees of freedom each (except for the neutrinos which in principle may only have two 
degrees of freedom) corresponding to a left or right-handed particle or anti-particle. 
The field labelled e contains a left-handed electron and a right-handed anti-electron, 
while e contains a left-handed anti-electron and a right-handed electron. Thus all four 
degrees of freedom are naturally (this is in accord with Lorentz invariance) contained in 
two independent fields e, e. This distinction is a property of Nature, since the charges 
of the SM particles depend on their handedness. In fact in each family we have five 
different charge multiplets given by 



where Q is a triplet under color SU(3), a doublet under weak SU(2) and carries U(l)y 
weak hypercharge Y — 1/3. The color anti-triplets (u, d) are singlets under SU{2) with 




(1) 
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Y — (— |, |) and finally the leptons appear as a electroweak doublet (L) and singlet (e) 
with Y = — 1,+2 respectively. Note, by definition, leptons are color singlets and thus 
do not feel the strong forces. The electric charge for all the quarks and leptons is given 
by the relation Qem = ^3 + \ where the (upper, lower) component of a weak doublet 
has T 3 = (+1/2, —1/2). Finally the Higgs boson multiplet, 

with Y = +1 is necessary to give mass to the W ± , Z° and to all quarks and leptons. In 
the SM vacuum, the field h° obtains a non-zero vacuum expectation value (h°) = v/y/2. 
Particle masses are then determined by the strength of the coupling to the Higgs. The 
peculiar values of the quark, lepton and Higgs charges is one of the central unsolved 
puzzles of the SM. The significance of this problem only becomes clear when one realizes 
that the interactions of all the particles (quarks, leptons, and Higgs bosons), via the 
strong and electroweak forces, are completely fixed by these charges. 

Let us now summarize the list of fundamental parameters needed to define the 
SM. If we do not include gravity or neutrino masses, then the SM has 19 fundamental 
parameters. These include the Z° and Higgs masses (Mz, va^) setting the scale for 
electroweak physics. The three gauge couplings Oii(Mz), % = 1,2,3, the 9 charged 
fermion masses and 4 quark mixing angles. Lastly, there is the QCD theta parameter 
which violates CP and thus is experimentally known to be less than 1CT 10 . Gravity 
adds one additional parameter, Newton's constant Gn = 1/M^ >1 or equivalently the 
Planck scale. Finally neutrino masses and mixing angles have been definitively observed 
in many recent experiments measuring solar and atmospheric neutrino oscillations, 
and by carefully measuring reactor or accelerator neutrino fluxes. The evidence for 
neutrino masses and flavor violation in the neutrino sector has little controversy. 
It is the first strong evidence for new physics beyond the SM. We shall return to 
these developments later. Neutrino masses and mixing angles are described by 9 new 
fundamental parameters - 3 masses, 3 real mixing angles and 3 CP violating phases. 

Let us now consider the minimal supersymmetric standard model [MSSM]. It 
is defined by the following two properties, (i) the particle spectrum , and (ii) their 
interactions. 

(i) Every matter fermion of the SM has a bosonic superpartner. In addition, every 
gauge boson has a fermionic superpartner. Finally, while the SM has one Higgs 
doublet, the MSSM has two Higgs doublets. 

= (3) 

with Y — +1, —1. The two Higgs doublets are necessary to give mass to up quarks, 
and to down quarks and charged leptons, respectively. The vacuum expectation 
values are now given by (h°) = vsin(3/\/2, (Ft ) = v cos (3/^/2 where tan (3 is a new 
free parameter of the MSSM. 



Desperately Seeking Supersymmetry [SUSY] 



6 



(ii) The MSSM has the discrete symmetry called R parity. || All SM particles are 
R parity even, while all superpartners are R odd. This has two important 
consequences. 

• The lightest superpartner [LSP] is absolutely stable, since the lightest state 
with odd R parity cannot decay into only even R parity states. Assuming that the 
LSP is electrically neutral, it is a weakly interacting massive particle. Hence it is a 
very good candidate for the dark matter of the universe. 

• Perhaps more importantly, the interactions of all superpartners with SM 
particles is completely determined by supersymmetry and the observed interactions 
of the SM. Hence, though we cannot predict the masses of the superpartners, we 
know exactly how they interact with SM particles. 

The MSSM has some very nice properties. It is perturbative and easily consistent 
with all precision electroweak data. In fact global fits of the SM and the MSSM provide 
equally good fits to the data [de Boer and Sander (2003)]. Moreover as the fiparticle 
masses increase, they decouple from low energy physics. On the other hand their masses 
cannot increase indefinitely since one soon runs into problems of "naturalness." In the 
SM the Higgs boson has a potential with a negative mass squared, of order the Z° mass, 
and an arbitrary quartic coupling. The quartic coupling stabilizes the vacuum value of 
the Higgs. In the MSSM the quartic coupling is fixed by supersymmetry in terms of the 
electroweak gauge couplings. As a result of this strong constraint, at tree level the light 
Higgs boson mass is constrained to be lighter than M z . One loop corrections to the 
Higgs mass are significant. Nevertheless the Higgs mass is bounded to be lighter than 
about 135 GeV [Okada et al (1991), Ellis et al (1991), Casas et al (1995), Carena et al 
(1995,1996), Haber et al (1997), Zhang (1999), Espinosa and Zhang (2000a,b), Degrassi 
et al (2003)]. The upper bound is obtained in the limit of large tan/3. 

It was shown early on that, even if the tree level Higgs mass squared was positive, 
radiative corrections due to a large top quark Yukawa coupling are sufficient to drive 
the Higgs mass squared negative [Ibanez and Ross (1982), Alvarez-Gaume et al (1983), 
Ibanez and Ross (1992)]. Thus radiative corrections naturally lead to electroweak 
symmetry breaking at a scale determined by squark and slepton Bbreaking masses. 
Note, a large top quark Yukawa coupling implies a heavy top quark. Early predictions 
for a top quark with mass above 50 GeV [Ibanez and Lopez (1983)] were soon challenged 
by the announcement of the discovery of the top quark by UAl with a mass of 40 GeV. 
Of course, this false discovery was much later followed by the discovery of the top quark 
at Fermilab with a mass of order 175 GeV. 

If the only virtue of Bis to explain why the weak scale (M^, rrih) is so much less 

|| One may give up R parity at the expense of introducing many new interactions with many new 
arbitrary couplings into the MSSM. These interactions violate either baryon or lepton number. Without 
R parity the LSP is no longer stable. There are many papers which give limits on these new couplings. 
The strongest constraint is on the product of couplings for the dimension four baryon and lepton 
number violating operators which contributes to proton decay. We do not discuss R parity violation 
further in this review. 
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than the Planck scale, one might ponder whether the benefits outweigh the burden of 
doubling the SM particle spectrum. Moreover there are many other ideas addressing 
the hierarchy problem, such as Technicolor theories with new strong interactions at a 
TeV scale. One particularly intriguing possibility is that the universe has more than 3 
spatial dimensions. In these theories the fundamental Planck scale M* is near a TeV, 
so there is no apparent hierarchy. I say apparent since in order to have the observed 
Newton's constant l/Mp t much smaller than 1/Mj 2 one needs a large extra dimension 
such that the gravitational lines of force can probe the extra dimension. If we live on a 3 
dimensional brane in this higher dimensional space then at large distances compared to 
the size of the d extra dimensions we will observe an effective Newton's constant given by 
G N = l/M% l = l/(R d M* +2 ) [Arkani-Hamed et al (1998)]. For example with d = 2 and 
M* = 1 TeV we need the radius of the extra dimension R ps 1 mm. If any of these new 
scenarios with new strong interactions at a TeV scaled are true then we should expect a 
plethora of new phenomena occurring at the next generation of high energy accelerators, 
i.e. the Large Hadron Collider [LHC] at CERN. It is thus important to realize that Bdoes 
much more. It provides a framework for understanding the 16 parameters of the SM 
associated with gauge and Yukawa interactions and also the 9 parameters in the neutrino 
sector. This will be discussed in the context of supersymmetric grand unified theories 
[SUSY GUTs ] and family symmetries. As we will see these theories are very predictive 
and will soon be tested at high energy accelerators or underground detectors. We will 
elaborate further on this below. Finally it is also naturally incorporated into string 
theory which provides a quantum mechanical description of gravity. Unfortunately this 
last virtue is apparently true for all the new ideas proposed to solve the gauge hierarchy 
problem. 

A possible subtitle for this article could be "A Tale of Two Symmetries: SUSY 
GUTs ." Whereas Bby itself provides a framework for solving the gauge hierarchy 
problem, i.e. why M z < M GUT , SUSY GUTs (with the emphasis on GUTs) adds 
the framework for understanding the relative strengths of the three gauge couplings and 
for understanding the puzzle of charge and mass. It also provides a theoretical lever 
arm for uncovering the physics at the Planck scale with experiments at the weak scale. 
Without any exaggeration it is safe to say that SUSY GUTs also address the following 
problems. 

• They explain charge quantization since weak hypercharge (Y) is imbedded in a 
non-abelian symmetry group. 

• They explain the family structure and in particular the peculiar color and 
electroweak charges of fermions in one family of quarks and leptons. 

• They predict gauge coupling unification. Thus given the experimentally determined 
values of two gauge couplings at the weak scale, one predicts the value of the third. 
The experimental test of this prediction is the one major success of Btheories. It 

Field theories in extra dimensions are divergent and require new non-perturbative physics, perhaps 
string theory, at the TeV scale. 
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relies on the assumption of fiparticles with mass in the 100 GeV to 1 TeV range. 
Hence it predicts the discovery of fiparticles at the LHC. 

• They predict Yukawa coupling unification for the third family. In SU(5) we obtain 
b — r unification, while in SO(10) we have t — b — r unification. We shall argue that 
the latter prediction is eminently testable at the Tevatron, the LHC or a possible 
Next Linear Collider. 

• With the addition of family symmetry they provide a predictive framework for 
understanding the hierarchy of fermion masses. 

• It provides a framework for describing the recent observations of neutrino masses 
and mixing. At zeroth order the See - Saw scale for generating light neutrino masses 
probes physics at the GUT scale. 

• The LSP is one of the best motivated candidates for dark matter. Moreover 
back of the envelope calculations of LSPs, with mass of order 100 GeV and 
annihilation cross-sections of order 1/TeV 2 , give the right order of magnitude 
of their cosmological abundance for LSPs to be dark matter. More detailed 
calculations agree. Underground dark matter detectors will soon probe the 
mass/cross-section region in the LSP parameter space. 

• Finally the cosmological asymmetry of baryons vs. anti-baryons can be explained 
via the process known as leptogenesis [Fukugita and Yanagida (1986)]. In this 
scenario an initial lepton number asymmetry, generated by the out of equilibrium 
decays of heavy Majorana neutrinos, leads to a net baryon number asymmetry 
today. 

Grand unified theories are the natural extension of the standard model. Ever 
since it became clear that quarks are the fundamental building blocks of all strongly 
interacting particles, protons, neutrons, pions, kaons, etc. and that they appear to be 
just as elementary as leptons, it was proposed [Pati and Salam (1973a,b, 1974)] that 
the strong SU(3) color group should be extended to SU(4) color with lepton number as 
the fourth color. 



In the Pati-Salam [PS] model, quarks and leptons of one family are united into two 
irreducible representations (Eqn. 0J). The two Higgs doublets of the MSSM sit in 
one irreducible representation (Eqn. EJ). This has significant consequences for fermion 
masses as we discuss later. However the gauge groups are not unified and there are 
still three independent gauge couplings, or two if one enlarges PS with a discrete parity 
symmetry where L R. PS must be broken spontaneously to the SM at some large 
scale Mq. Below the PS breaking scale the three low energy couplings c^, i = 1,2,3 



Gp a ti — Salam — 

SU 4 (color) x SU 2 (L) x SU 2 (R) 




(4) 



(H u H d ) 



(5) 
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renormalize independently. Thus with the two gauge couplings and the scale Mq one 
can fit the three low energy couplings. 

Shortly after PS, the completely unified gauge symmetry SU5 was proposed 
[Georgi and Glashow (1974)]. Quarks and leptons of one family sit in two irreducible 
represent at ions . 



{Q = " e c u c } C 10, (6) 





K L ={ e p c 5 ( 7 ) 

The two Higgs doublets necessarily receive color triplet SU(5) partners filling out 5, 5 
representations. 

C 5 H , 5 h (8) 

As a consequence of complete unification the three low energy gauge couplings are given 
in terms of only two independent parameters, the one unified gauge coupling o>g(Mc) 
and the unification (or equivalently the SU(5) symmetry breaking ) scale Mq [Georgi 
et al (1974)]. Hence there is one prediction. In addition we now have the dramatic 
prediction that a proton is unstable to decay into a n° and a positron, e + . 

Finally complete gauge and family unification occurs in the group SO10 [Georgi 
(1974), Fritzsch and Minkowski (1974)] with one family contained in one irreducible 
representation 

10 + 5 + u C 16 (9) 
and the two multiplets of Higgs unified as well. 

5 H , 5 h C 10 h , (10) 

(See Tabled). 

GUTs predict that protons decay with a lifetime t p of order M G /a G m? p . The first 
experiments looking for proton decay were begun in the early 1980s. However at the 
very moment that proton decay searches began, motivated by GUTs, it was shown that 
SUSY GUTs naturally increase Mq, thus increasing the proton lifetime. Hence, if SUSY 
GUTs were correct, it was unlikely that the early searches would succeed [Dimopoulos 
et al (1981), Dimopoulos and Georgi (1981), Ibanez and Ross (1981), Sakai (1981), 
Einhorn and Jones (1982), Marciano and Senjanovic (1982)]. At the same time, it was 
shown that SUSY GUTs did not significantly affect the predictions for gauge coupling 
unification (for a review see [Dimopoulos et al (1991), Raby (2002a)]). At present, 
non- SUSY GUTs are excluded by the data for gauge coupling unification; where as 
SUSY GUTs work quite well. So well in fact, that the low energy data is now probing 
the physics at the GUT scale. In addition, the experimental bounds on proton decay 
from Super-Kamiokande exclude non-BGUTs, while severely testing BGUTs. Moreover, 



Desperately Seeking Supersymmetry [SUSY] 



10 



Table 1. This table gives the particle spectrum for the 16 dimensional spinor 
representation of SO(10). The states are described in terms of the tensor product 
of five spin 1/2 states with spin up (+) or down (— ) and in addition having an even 
number of (— ) spins. { r,b,y } are the three color quantum numbers, and Y is weak 
hypercharge given in terms of the formula §£(C) — S(W) where the sum (£) is over all 
color and weak spins with values (± |). Note, an SO(10) rotation corresponds either 
to raising one spin and lowering another or raising (or lowering) two spins. In the table, 
the states are arranged in SU(5) multiplets. One readily sees that the first operation of 
raising one spin and lowering another is an SU(5) rotation, while the others are special 
to SO(10). 



State Y Color Weak 

= |S(C) - S(W) C spins W spins 
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future underground proton decay/neutrino observatories, such as the proposed Hyper- 
Kamiokande detector in Japan or UNO in the USA will cover the entire allowed range 
for the proton decay rate in SUSY GUTs . 

If Bis so great, if it is Nature, then where are the Bparticles? Experimentalists at 
high energy accelerators, such as the Fermilab Tevatron and the CERN LHC (now under 
construction), are desperately seeking Bparticles or other signs of B. At underground 
proton decay laboratories, such as Super-Kamiokande in Japan or Soudan II in 
Minnesota, USA, electronic eyes continue to look for the tell-tale signature of a proton 
or neutron decay. Finally, they are searching for cold dark matter, via direct detection 
in underground experiments such as CDMS, UKDMC or EDELWEISS, or indirectly by 
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searching for energetic gammas or neutrinos released when two neutralino dark matter 
particles annihilate. In Sect. 121 we focus on the perplexing experimental/theoretical 
problem of where are these Bparticles. We then consider the status of gauge coupling 
unification (Sect. OJ), proton decay predictions (Sect. 0J, fermion masses and mixing 
angles (including neutrinos) and the fiflavor problem (Sect. EJ), and fidark matter (Sect. 
We conclude with a discussion of some remaining open questions. 

2. Where are the supersymmetric particles ? 

The answer to this question depends on two interconnected theoretical issues - 

(i) the mechanism for fibreaking, and 

(ii) the scale of fibreaking. 

The first issue is inextricably tied to the Bflavor problem. While the second issue is tied 
to the gauge hierarchy problem. We discuss these issues in sections 12.11 and 12.21 

2.1. fiBreaking Mechanisms 

Supersymmetry is necessarily a local gauge symmetry, since Einstein's general theory of 
relativity corresponds to local Poincare symmetry and supersymmetry is an extension 
of the Poincare group. Hence Bbreaking must necessarily be spontaneous, in order not 
to cause problems with unitarity and/or relativity. In this section we discuss some 
of the spontaneous Bbreaking mechanisms considered in the literature. However from 
a phenomenological stand point, any spontaneous fibreaking mechanism results in soft 
fibreaking operators with dimension 3 or less (such as quadratic or cubic scalar operators 
or fermion mass terms) in the effective low energy theory below the scale of fibreaking 
[Dimopoulos and Georgi (1981), Sakai (1981), Girardello and Grisaru (1982)]. There 
are a priori hundreds of arbitrary soft fibreaking parameters (the coefficients of the 
soft fibreaking operators) [Dimopoulos and Sutter (1995)]. These are parameters not 
included in the SM but are necessary to compare with data or make predictions for new 
experiments. 

The general set of renormalizable soft fibreaking operators, preserving the solution 
to the gauge hierarchy problem, is given in a paper by [Girardello and Grisaru (1982)]. 
These operators are assumed to be the low energy consequence of spontaneous fibreaking 
at some fundamental fibreaking scale Ag ^> Mz- The list of soft fibreaking parameters 
includes squark and slepton mass matrices, cubic scalar interaction couplings, gaugino 
masses, etc. Let us count the number of arbitrary parameters [Dimopoulos and Sutter 
(1995)]. Left and right chiral scalar quark and lepton mass matrices are a priori 
independent 3x3 hermitian matrices. Each contains 9 arbitrary parameters. Thus 
for the scalar partners of {Q, u, d, L, e} we have 5 such matrices or 45 arbitrary 
parameters. In addition corresponding to each complex 3x3 Yukawa matrix (one for 
up and down quarks and charged leptons) we have a complex soft fibreaking trilinear 
scalar coupling (A) of left and right chiral squarks or sleptons to Higgs doublets. This 
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adds 3 x 18 = 54 additional arbitrary parameters. Finally, add to these 3 complex 
gaugino masses (M i? % = 1,2,3), and the complex soft Bbreaking scalar Higgs mass 
(fiB) and we have a total of 107 arbitrary soft Bbreaking parameters. In additon, the 
minimal supersymmetric extension of the SM requires a complex Higgs mass parameter 
(//) which is the coefficient of a supersymmetric term in the Lagrangian. Therefore, 
altogether this minimal extension has 109 arbitrary parameters. Granted, not all of 
these parameters are physical. Just as not all 54 parameters in the three complex 
3x3 Yukawa matrices for up and down quarks and charged leptons are observable. 
Some of them can be rotated away by unitary redefinitions of quark and lepton 
superfields. Consider the maximal symmetry of the kinetic term of the theory — global 
SU(3) Q x SU(3) a x SU(3)a x SU(3) L x SU(3) e x U(lf x U(1) R . Out of the total 
number of parameters - 163 = 109 (new Bparameters) + 54 (SM parameters) - we can 
use the SU(3) 5 to eliminate 40 parameters and 3 of the U(l)s to remove 3 phases. The 
other 3 U(l)s however, are symmetries of the theory corresponding to B, L and weak 
hypercharge Y. We are thus left with 120 observables, corresponding to the 9 charged 
fermion masses, 4 quark mixing angles and 107 new, arbitrary observable fiparameters. 

Such a theory, with so many arbitrary parameters, clearly makes no predictions. 
However, this general MSSM is a "straw man" (one to be struck down), but fear not 
since it is the worst case scenario. In fact, there are several reasons why this worst case 
scenario cannot be correct. First, and foremost, it is severely constrained by precision 
electroweak data. Arbitrary 3x3 matrices for squark and slepton masses or for trilinear 
scalar interactions maximally violate quark and lepton flavor. The strong constraints 
from flavor violation were discussed by [Dimopoulos and Georgi (1981), Dimopoulos 
and Sutter (1995), Gabbiani et al (1996)]. In general, they would exceed the strong 
experimental contraints on flavor violating processes, such as b — > sj, or b — > sl + l~, 
B s — > ji — > ej, /i — e conversion in nuclei, etc. In order for this general MSSM 

not to be excluded by flavor violating constraints, the soft Bbreaking terms must be 
either 

(i) flavor independent, 

(ii) aligned with quark and lepton masses or 

(iii) the first and second generation squark and slepton masses (m) should be large (i.e. 
greater than a TeV). 

In the first case, squark and slepton mass matrices are proportional to the 3x3 identity 
matrix and the trilinear couplings are proportional to the Yukawa matrices. In this 
case the squark and slepton masses and trilinear couplings are diagonalized in the same 
basis that quark and lepton Yukawa matrices are diagonalized. This limit preserves 
three lepton numbers - L e , L^, L T - (neglecting neutrino masses) and gives minimal 
flavor violation (due only to CKM mixing) in the quark sector [Hall et al (1986)]. 
The second case does not require degenerate quark flavors, but approximately diagonal 
squark and slepton masses and interactions, when in the diagonal quark and lepton 
Yukawa basis. It necessarily ties any theoretical mechanism explaining the hierarchy 
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of fermion masses and mixing to the hierarchy of sfermion masses and mixing. This 
will be discussed further in Sections 15.21 and 15.41 Finally, the third case minimizes 
flavor violating processes, since all such effects are given by effective higher dimension 
operators which scale as 1/m 2 . The theoretical issue is what Bbreaking mechanisms are 
"naturally" consistent with these conditions. 

Several such Bbreaking mechanisms exist in the literature. They are called minimal 
supergravity [mSugra] breaking, gauge mediated Bbreaking [GMSB], dilaton mediated 
[DMSB], anomaly mediated [AMSB], and gaugino mediated [gMSB]. Consider first 
mSugra which has been the benchmark for experimental searches. The minimal 
supergravity model [Ovrut and Wess (1982), Chamseddine et al (1982), Barbieri et al 
(1982), Ibanez (1982), Nilles et al (1983), Hall et al (1983)] is defined to have the minimal 
set of soft Bbreaking parameters. It is motivated by the simplest ([Polony (1977)] 
hidden sector in supergravity with the additional assumption of grand unification. This 
Bbreaking scenario is also known as the constrained MSSM [CMSSM] [Kane et al (1994)]. 
In mSUGRA/CMSSM there are four soft Bbreaking parameters at Mq defined by mo, 
a universal scalar mass; A , a universal trilinear scalar coupling; M1/2, a universal 
gaugino mass; and /io-B , the s °ft Bbreaking Higgs mass parameter where /zo, is the 
super symmetric Higgs mass parameter. In most analyses, | /x 1 an d HqBq are replaced, 
using the minimization conditions of the Higgs potential, by Mz and the ratio of Higgs 
VEVs tan/3 = (H u )/(Hd). Thus the parameter set defining mSugra/CMSSM is given 
by 

m , A , Mi/2, sign(/i ), tan/3. (11) 

This scenario is an example of the first case above (with minimal flavor violation), 
however it is certainly not a consequence of the most general supergravity theory and 
thus requires further justification. Nevertheless it is a useful framework for experimental 
Bsearches. 

In GMSB models, Bbreaking is first felt by messengers carrying standard model 
charges and then transmitted to to the superpartners of SM particles [sparticles] via 
loop corrections containing SM gauge interactions. Squark and slepton masses in these 
models are proportional to a« with = F/M. In this expression, cti, i = 1,2,3 
are the fine structure constants for the SM gauge interactions, y/~F is the fundamental 
scale of Bbreaking, M is the messenger mass, and Ag is the effective Bbreaking scale. 
In GMSB the flavor problem is naturally solved since all squarks and sleptons with the 
same SM charges are degenerate and the A terms vanish to zeroth order. In addition, 
GMSB resolves the formidable problems of model building [Fayet and Ferrara (1977)] 
resulting from the direct tree level Bbreaking of sparticles. This problem derives from 
the supertrace theorem, valid for tree level Bbreaking, 

S(2J+ 1)(-1) 2J M) = (12) 

where the sum is over all particles with spin J and mass Mj. It generically leads 
to charged scalars with negative mass squared [Fayet and Ferrara (1977), Dimopoulos 
and Georgi (1981)]. Fortunately the supertrace theorem is explicitly violated when 
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Bbreaking is transmitted radiatively. + Low energy Bbreaking models [Dimopoulos and 
Raby (1981), Dine et al (1981), Witten (1981), Dine and Fischler (1982), Alvarez- 
Gaume et al (1982)], with \[F ~ M ~ A 6 ~ 100 TeV make dramatic predictions 
[Dimopoulos et al (1996)]. Following the seminal work of [Dine and Nelson (1993), Dine 
et al (1995,1996)] complete GMSB models now exist (for a review, see [Giudice and 
Rattazzi (1999)]. Of course the fundamental Bbreaking scale can be much larger than 
the weak scale. Note Bbreaking effects are proportional to 1/M and hence they decouple 
as M increases. This is a consequence of Bbreaking decoupling theorems [Polchinski 
and Susskind (1982), Dimopoulos and Raby (1983), Banks and Kaplunovsky (1983)]. 
However when \J~F~ > 10 10 GeV then supergravity becomes important. 

DMSB, motivated by string theory, and AMSB and gMSB, motivated by brane 
models with extra dimensions, also alleviate the Bflavor problem. We see that 
there are several possible Bbreaking mechanisms which solve the Bflavor problem and 
provide predictions for superpartner masses in terms of a few fundamental parameters. 
Unfortunately we do not a priori know which one of these (or some other) Bbreaking 
mechanism is chosen by nature. For this we need experiment. 

2.2. Fine Tuning or "Naturalness" 

Presently, the only evidence for supersymmetry is indirect, given by the successful 
prediction for gauge coupling unification. Supersymmetric particles at the weak scale are 
necessary for this to work, however it is discouraging that there is yet no direct evidence. 
Searches for new supersymmetric particles at CERN or Fermilab have produced only 
lower bounds on their mass. The SM Higgs mass bound, applicable to the MSSM when 
the CP odd Higgs (A) is much heavier, is 114.4 GeV [LEP2 (2003)]. In the case of an 
equally light A, the Higgs bound is somewhat lower ~ 89 GeV. Squark, slepton and 
gluino mass bounds are of order 200 GeV, while the chargino bound is 103 GeV [LEP2 
(2003)]. In addition other indirect indications for new physics beyond the standard 
model, such as the anomalous magnetic moment of the muon (a M ), are inconclusive. 
Perhaps Nature does not make use of this beautiful symmetry? Or perhaps the Bparticles 
are heavier than we once believed. 

Nevertheless, global fits to precision electroweak data in the SM or in the MSSM 
give equally good x 2 /d°f [de Boer (2003)]. In fact the fit is slightly better for the 
MSSM due mostly to the pull of a M . The real issue among Benthusiasts is the problem 
of fine tuning. If Bis a solution to the gauge hierarchy problem (making the ratio 
Mz/Mq ~ 10~ 14 "naturally" small), then radiative corrections to the Z mass should 
be insensitive to physics at the GUT scale, i.e. it should not require any "unnatural" 
fine tuning of GUT scale parameters. A numerical test of fine tuning is obtained by 
defining the fine tuning parameter A a (= \jfr-^-\), the logarithmic derivative of the Z 
mass with respect to different "fundamental" parameters a a = {/i, Mi/ 2 , mo, A , . . .} 

+ It is also violated in supergravity where the right-hand side is replaced by 2(N — l)m^ 2 with m 3 / 2} 
the gravitino mass and N the number of chiral superfields. 
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defined at M G [Ellis et al (1986), Barbieri and Giudice (1988), de Carlos and Casas 
(1993), Anderson et al (1995)]. Smaller values of A a correspond to less fine tuning and 
roughly speaking p = Max(A Q ,)~ 1 is the probability that a given model is obtained in 
a random search over Bparameter space. 

There are several recent analyses, including LEP2 data, by [Chankowski et al 
(1997), Barbieri and Strumia (1998), Chankowski et al (1999)]. In particular [Barbieri 
and Strumia (1998), Chankowski et al (1999)] find several notable results. In their 
analysis [Barbieri and Strumia (1998)] only consider values of tan/? < 10 and soft 
Bbreaking parameters of the CMSSM or gauge-mediated Bbreaking. [Chankowski et al 
(1999)] also consider large tan/5 = 50 and more general soft Bbreaking scenarios. They 
both conclude that the value of Max(A a ) is significantly lower when one includes the 
one loop radiative corrections to the Higgs potential as compared to the tree level Higgs 
potential used in the analysis of [Chankowski et al (1997)]. In addition they find that the 
experimental bound on the Higgs mass is a very strong constraint on fine tuning. Larger 
values of the light Higgs mass require larger values of tan/5. Values of Max(A a ) < 10 
are possible for a Higgs mass < 111 GeV (for values of tan (3 < 10 used in the analysis of 
[Barbieri and Strumia (1998)]). However allowing for larger values of tan (3 [Chankowski 
et al (1999)] allows for a heavier Higgs. With LEP2 bounds on a SM Higgs mass of 114.4 
GeV, larger values of tan (3 > 4 are required. It is difficult to conclude too much from 
these results. Note, the amount of fine tuning is somewhat sensitive to small changes 
in the definition of A a . For example, replacing a a — > a 2 a or M| — > Mz will change the 
result by a factor of 2 ±1 . Hence factors of two in the result are definition dependent. 
Let us assume that fine tuning by 1/10 is acceptable, then is fine tuning by 1 part in 
100 "unnatural." Considering the fact that the fine tuning necessary to maintain the 
gauge hierarchy in the SM is at least 1 part in 10 28 , a fine tuning of 1 part in 100 (or 
even 10 3 ) seems like a great success. 

A slightly different way of addressing the fine tuning question says if I assign equal 
weight to all "fundamental" parameters at Mq and scan over all values within some 
a priori assigned domain, what fraction of this domain is already excluded by the low 
energy data. This is the analysis that [Strumia (1999)] uses to argue that 95% of the 
Bparameter space is now excluded by LEP2 bounds on the Bspectrum and in particular 
by the Higgs and chargino mass bounds. This conclusion is practically insensitive to 
the method of Bbreaking assumed in the analysis which included the CMSSM, gauge- 
mediated or anomaly mediated Bbreaking or some variations of these. One might still 
question whether the a priori domain of input parameters, upon which this analysis 
stands, is reasonable. Perhaps if we doubled the input parameter domain we could find 
acceptable solutions in 50% of parameter space. To discuss this issue in more detail, let 
us consider two of the parameter domains considered in [Strumia (1999)]. Within the 
context of the CMSSM, he considers the domain defined by 

m , |A)|, |Afi /2 |, |/i |, | -B 1 = (0-r 1) m fi (13) 
where (a -j- b) = a random number between a and b and the overall mass scale mu is 
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obtained from the minimization condition for electroweak symmetry breaking. He also 
considered an alternative domain defined by 

m = (-v3)m fil \fi \, |M 1/2 | = (- 3) m fi , (14) 

A , B = (-3 3) m 

with the sampling of mo, My 2 , /io using a flat distribution on a log scale. In both 
cases, he concludes that 95% of parameter space is excluded with the light Higgs and 
chargino mass providing the two most stringent constraints. Hence we have failed to 
find Bin 95% of the allowed region of parameter space. But perhaps we should open the 
analysis to other, much larger regions, of fiparameter space. We return to this issue in 
Sections O and O 

For now however, let us summarize our discussion of naturalness constraints with 
the following quote from [Chankowski et al (1999)], "We re-emphasize that naturalness 
is [a] subjective criterion, based on physical intuition rather than mathematical rigour. 
Nevertheless, it may serve as an important guideline that offers some discrimination 
between different theoretical models and assumptions. As such, it may indicate which 
domains of parameter space are to be preferred. However, one should be very careful 
in using it to set any absolute upper bounds on the spectrum. We think it safer to 
use relative naturalness to compare different scenarios, as we have done in this paper." 
As these authors discuss in their paper, in some cases the amount of fine tuning can 
be dramatically decreased if one assumes some linear relations between GUT scale 
parameters. These relations may be due to some, yet unknown, theoretical relations 
coming from the fundamental physics of Bbreaking, such as string theory. 

In the following we consider two deviations from the simplest definitions of fine 
tuning and naturalness. The first example, called focus point [FP] [Feng and Moroi 
(2000), Feng et al (2000a,b,c), Feng and Matchev (2001)] is motivated by infra-red fixed 
points of the renormalization group equations for the Higgs mass and other dimensionful 
parameters. The second had two independent motivations. In the first case it is 
motivated by the fiflavor problem and in this incarnation it is called the radiative 
inverted scalar mass hierarchy [RISMH] [Bagger et al (1999,2000)]. More recently it 
was reincarnated in the context of SO (10) Yukawa unification for the third generation 
quarks and leptons [YU] [Raby (2001), Dermfsek (2001), Baer and Ferrandis (2001), 
Blazek et al (2002a,b), Auto et al (2003), Tobe and Wells (2003)] scenarios. In both 
scenarios the upper limit on soft scalar masses is increased much above 1 TeV. 

2. 3. The focus point region of fibreaking parameter space 

In the focus point Bbreaking scenario, Feng et al [Feng and Moroi (2000), Feng et 
al (2000a,b)] consider the renormalization group equations [RGE] for soft Bbreaking 
parameters, assuming a universal scalar mass mo at Mq- This may be as in the CMSSM 
(Eqn. n"Tj) or even a variation of AMSB. They show that, if the top quark mass is 
approximately 174 GeV, then the RGEs lead to a Higgs mass which is naturally of 
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order the weak scale, independent of the precise value of m which could be as large as 
3 TeV. It was also noted that the only fine tuning in this scenario was that necessary to 
obtain the top quark mass, i.e. if the top quark mass is determined by other physics then 
there is no additional fine tuning needed to obtain electroweak symmetry breaking.* As 
discussed in [Feng et al (2000c)] this scenario opens up a new window for neutralino 
dark matter. Cosmologically acceptable neutralino abundances are obtained even with 
very large scalar masses. Moreover as discussed in [Feng and Matchev (2001)] the focus 
point scenario has many virtues. In the limit of large scalar masses, gauge coupling 
unification requires smaller threshold corrections at the GUT scale, in order to agree 
with low energy data. In addition, larger scalar masses ameliorate the Bflavor and 
CP problems. This is because both processes result from effective higher dimensional 
operators suppressed by two powers of squark and/or slepton masses. Finally a light 
Higgs mass in the narrow range from about 114 to ~ 120 GeV is predicted. Clearly the 
focus point region includes a much larger range of soft Bbreaking parameter space than 
considered previously. It may also be perfectly "natural." 

The analysis of the focus point scenario was made within the context of the CMSSM. 
The focus point region extends to values of m up to 3 TeV. This upper bound increases 
from 3 to about 4 TeV as the top quark mass is varied from 174 to 179 GeV. On the 
other hand, as tan f3 increases from 10 to 50, the allowed range in the ijiq — Mi/ 2 plane 
for Aq = 0, consistent with electroweak symmetry breaking, shrinks. As we shall see 
from the following discussion, this narrowing of the focus point region is most likely an 
artifact of the precise CMSSM boundary conditions used in the analysis. In fact the 
CMSSM parameter space is particularly constraining in the large tan/3 limit. 

2.4- SO(10) Yukawa unification and the radiative inverted scalar mass hierarchy 
[RISMHJ 

The top quark mass M t ~ 174 GeV requires a Yukawa coupling X t ~ 1. In the minimal 
SO(10) fimodel [MSOi SM] the two Higgs doublets, H u , H d , of the MSSM are contained 
in one 10. In addition the three families of quarks and leptons are in 16j, i — 1, 2, 3. In 
the MSOioSM the third generation Yukawa coupling is given by 

A 16 3 10 16 3 = A (Q 3 H u t + L 3 H u v T + Q 3 H d b + L 3 H d f ) (15) 

+ A (^Qs Q 3 + t f) T + X(Q 3 L 3 + tb) T. 

Thus we obtain the unification of all third generation Yukawa couplings with 

A = X t = \ b = \ T = A„ T . (16) 

Of course this simple Yukawa interaction, with the constant A replaced by a 3 x 3 Yukawa 
matrix, does not work for all three families. Jj In this discussion, we shall assume that 

* For a counter discussion of fine tuning in the focus point region, see [Romanino and Strumia (2000)]. 
jj In such a theory there is no CKM mixing matrix and the down quark and charged lepton masses 
satisfy the bad prediction ^ ^ — = — ~ 2M • 
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the first and second generations obtain mass using the same 10, but via effective higher 
dimensional operators resulting in a hierarchy of fermion masses and mixing. In this 
case, Yukawa unification for the third family (Eqn. ITH|) is a very good approximation. 
The question then arises, is this symmetry relation consistent with low energy data 
given by 

M t = 174.3 ±5.1 GeV, m b (m b ) = 4.20 ± 0.20 GeV, (17) 
M T = 1.7770 ± 0.0018, 

where the error on M T is purely a theoretical uncertainty due to numerical errors in 
the analysis. Although this topic has been around for a long time, it is only recently 
that the analysis has included the complete one loop threshold corrections at the weak 
scale [Raby (2001), Derrmsek (2001), Baer and Ferrandis (2001), Blazek et al (2002a,b), 
Auto et al (2003), Tobe and Wells (2003)]. It turns out that these corrections are very 
important. The corrections to the bottom mass are functions of squark and gaugino 
masses times a factor of tan/? ~ m t (m t ) /mb(m t ) ~ 50. For typical values of the 
parameters the relative change in the bottom mass dm^/mb is very large, of order 
50%. At the same time, the corrections to the top and tau masses are small. For 
the top, the same one loop corrections are proportional to 1/ tan/3, while for the tau, 
the dominant contribution from neutralino loops is small. These one loop radiative 
corrections are determined, through their dependence on squark and gaugino masses, 
by the soft Bbreaking parameters at Ma- In the MSOioSM we assume the following 
dimensionful parameters. 

mie, mio, A»h, A), M 1/2 , /i (18) 

where mi 6 is the universal squark and slepton mass; 7?7# u/d = m io(lT^ m fl") is the Higgs 
up/down mass squared; A is the universal trilinear Higgs - scalar coupling; M 1 / 2 is the 
universal gaugino mass and \i is the super symmetric Higgs mass, tan (3 ps 50 is fixed by 
the top, bottom and r mass. Note, there are two more parameters (mio, Am H ) than 
in the CMSSM. They are needed in order to obtain electroweak symmetry breaking 
solutions in the region of parameter space with uiiq 3> /i, M1/2. We shall defer a more 
detailed discussion of the results of the MSOi SM to Sects. 15.11 and IdI Suffice it to say 
here that good fits to the data are only obtained in a narrow region of fiparameter space 
given by 

Aq PS —2 777.16, m 10 ~ 777i6, Amjj ^ 0.13 (19) 

7«i6 > /i ~ M1/2, mie > 2 TeV. 

Once more we are concerned about fine-tuning with mi6 so large. However, we 
discover a fortuitous coincidence. This region of parameter space (Eqn. ITHjl "naturally" 
results in a radiative inverted scalar mass hierarchy with fh\ 2 3> fh\ [Bagger et al 
(1999,2000)], i.e. first and second generation squark and slepton masses are of order m^ 6 , 
while the third generation scalar masses are much lighter. Since the third generation 
has the largest couplings to the Higgs bosons, they give the largest radiative corrections 
to the Higgs mass. Hence with lighter third generation squarks and sleptons, a light 
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Higgs is more "natural." Although a detailed analysis of fine-tuning parameters is not 
available in this regime of parameter space, the results of several papers suggest that 
the fine-tuning concern is minimal (see for example, [Dimopoulos and Giudice (1995), 
Chankowski et al (1999), Kane et al (2003)]). While there may not be any fine-tuning 
necessary in the MSOioSM region of Bparameter space (Eqn. there is still one 

open problem. There is no known fibreaking mechanism which "naturally" satisfies the 
conditions of Eqn. EH On the other hand, we conclude this section by noting that the 
latter two examples suggest that there is a significant region of Bbreaking parameter 
space which is yet to be explored experimentally. 



3. Gauge coupling unification 

The apparent unification of the three gauge couplings at a scale of order 3 x 10 16 GeV is, 
at the moment, the only experimental evidence for low energy supersymmetry [Amaldi 
et al (1991), Ellis et al (1991), Langacker and Luo (1991)]. In this section we consider 
the status of gauge coupling unification and the demise of minimal BSU(5). 

The theoretical analysis of unification is now at the level requiring two loop 
renormalization group running from Mq to Mz- Consistency then requires including 
one loop threshold corrections at both the GUT and weak scales. Once GUT threshold 
corrections are considered, a precise definition of the GUT scale (Mq) is needed. The 
three gauge couplings no longer meet at one scale, f since 

aiifj.)- 1 = a G 1 + A i ( f i) (20) 

where the corrections Aj(/i) are logarithmic functions of mass for all states with GUT 
scale mass. In principle, the GUT scale can now be defined as the mass Mj of the X, X 
gauge bosons mediating proton decay or as the scale where any two couplings meet. We 
define Mq as the value of \x where 

a x (M G ) = a 2 (M G ) = a G . (21) 
Using two loop RGE from Mz to Mq, we find 

Mq ~ 3 x 10 16 GeV (22) 
cyq 1 « 24. 

In addition, good fits to the low energy data require 

e 3 = MMg) ~ &g) ~ -3% to - 4%. (23) 

OLQ 

Note the exact value of the threshold correction (e 3 ), needed to fit the data, depends 
on the weak scale threshold corrections and in particular on the fiparticle spectrum. 
We shall return to this later. On the other hand, significant contributions to the GUT 
threshold correction €3 typically arise from the Higgs and GUT breaking sectors of the 

f [Brodsky et al (2003)] has argued that the three gauge couplings always meet in a GUT at a scale 
above the largest GUT mass. He defines this to be the GUT scale. Unfortunately, this scale cannot be 
defined in the effective low energy theory. 
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theory. Above Mq there is a single coupling constant olq ~ «g which then runs up to 
some fundamental scale M*, such as the string scale, where the running is cut off. The 
GUT symmetry, in concert with supersymmetry, regulates the radiative corrections. 
Without the GUT, e 3 would naturally take on a value of order one. 

Following [Lucas and Raby (1996)] we show that the allowed functional dependence 
of e 3 on GUT symmetry breaking vacuum expectation values [VEVs] is quite restricted. 
Consider a general SO(10) theory with 

(ni 6 + 3) 16 + me 16 + n 10 10 + n 45 45 + n 54 54 + n x 1. (24) 

Note, the superpotential for the GUT breaking sector of the theory typically has a U(l) n 
x R symmetry which, as we shall see, has an important consequence for the threshold 
corrections. The one loop threshold corrections are given by 

a; 1 (M G )=a G 1 -A l (25) 

with 
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The sum is over all super heavy particles £ with mass M^ and b\ is the contribution the 
super heavy particle would make to the beta function coefficient 6, if the particle were 
not integrated out at Mq. 

As a consequence of Band the U(l) symmetries, Lucas and Raby proved the 
following theorem: e 3 is only a function of U(l) and R invariant products of powers 
of VEVs {&}, i.e. 



6 3 = F(Ci,...,C m ) + ^log 



M'rp 



iff 



+ ■■■. (27) 



As an example, consider the symmetry breaking sector given by the superpotential 

W symbreaking = T^^iC^l + ^S^l + ^Ai^) (28) 

+ A 2 (^ + s 1 A) + sA 2 

+ S'(SS 2 + A 1 A)+S 3 S' 2 

where the fields transform as follows {A 1: A 2 , A, A[} C 45, {S, S'} C 54, 

t/>, V> C 16, 16, and {iSi, • • • , 1S4} C 1. In addition we include the Lagrangian for 

the electroweak Higgs sector given by 

Lmggs = [IO1A1IO2 + S 5 10 2 2 \ F + ^[z*10j\ D . (29) 

Wsym breaking has a [U(l)] 4 x R symmetry. Since Bis unbroken, the potential has 
many F and D flat directions. One in particular (Eqn. l3*Uj) breaks SO(10) to 
SU(3)xSU(2)xU(l)y leaving only the states of the MSSM massless plus some non- 
essential SM singlets. 

(A 1 ) = a 1 ^(B-L), (A 2 ) = a 2 h (30) 
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Table 2. Recent preliminary lower bounds on the proton lifetime into specific decay 
modes from Super-Kamiokande [Jung (2002)]. 

mode r/B limit [10 33 yrs] 

p^n a + e+ 5.7 @ 90%C.L. 
p -> K+ + v 1.9 @ 90%C.L. 



2x2 



(A) (S) = s diag(lll - 3/2 - 3/2) ® 1 

(</,) = u |SU(5) singlet) ($) = U |SU(5) singlet) 

The VEVs {ax, a, v,v, S4} form a complete set of independent variables along the F 
and D flat directions. Note since the superpotential (Eqn. |2*H|) contains higher dimension 
operators fixed by the cutoff scale M* ~ Mpi anc k the GUT scale spectrum ranges from 
10 13 — 10 20 GeV. Nevertheless the threshold corrections are controlled. The only invariant 

a 4 

under a [U(l)j xR rotation of the VEVs is ( = ^a^r. By an explicit calculation we find 
the threshold correction 



3a G r 256 1 , 
= < log log 

d 5tt I 5 243 2 5 



(1 - 25C)' 



(i-C) 

Taking reasonable values of the VEVs given by 

CL\ = 2a,2 — 21S4 = Mq 
and the effective color triplet Higgs mass 



+ log 



M e T ff 



M, 



G 



M% n ~ 10 19 GeV 



we find 



(31) 



(32) 



(33) 



(34) 



£ = 16 e 3 w -0.030. 

Note, the large value of My is necessary to suppress proton decay rates as discussed 
in the following section. 



4. Nucleon Decay : Minimal SU(5) fiGUT 

Protons and neutrons [nucleons] are not stable particles; they necessarily decay in any 
GUT. Super-Kamiokande and Soudan II are looking for these decay products. The most 
recent (preliminary) Super-Kamiokande bounds on the proton lifetime [Jung (2002)] 
are given in Table [21 In the future, new detectors > 10 times larger than Super-K 
have been proposed - Hyper-Kamiokande in Japan and UNO in the USA. Note, a 
generic, dimension 6 nucleon decay operator is given by a 4 Fermion operator of the 
form ~ (1/A 2 ) q q q I. Given the bound r p > 5 x 10 33 yrs we find A > 4 x 10 15 GeV. 
This is nice, since it is roughly consistent with the GUT scale and with the See-Saw 
scale for neutrino masses. 

In this section we consider nucleon decay in the Minimal SU(5) BGUT in more 
detail. In minimal BSU(5), we have the following gauge and Higgs sectors. The gauge 
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sector includes the gauge bosons for SU(5) which decompose, in the SM, to SU(3) x 
SU(2) x U(l) plus the massive gauge bosons {X, X}. The {X, X} bosons with charges 
{(3,2, —5/3), (3,2,5/3)} are responsible for nucleon decay. The minimal SU(5) theory- 
has, by definition, the minimal Higgs sector. It includes a single adjoint of SU(5), E C 
24, for the GUT breaking sector and the electroweak Higgs sector (Eqn. |HJ) 

C 5 H , 5 H - 

The superpotential for the GUT breaking and Higgs sectors of the model is given by 
[Witten (1981), Dimopoulos and Georgi (1981), Sakai (1981)] 

W = - Tr S 3 + — Tr S 2 + H(E + W)H. (35) 
3 2 

In general, nucleon decay can have contributions from operators with dimensions 4, 5 
and 6. 




4-1. Dimension 6 operators 

The dimension 6 operators are derived from gauge boson exchange (see Fig. [T}. We 
obtain the effective dimension 6 (four Fermion) operator given by 

Jf # « s L - ' 36 > 

Thus the decay amplitude is suppressed by one power of \jM\. How is M x related 
to the GUT scale M G determined by gauge coupling unification? Recall, in general we 
have 

e 3 = ^ In ^ + 6 3 (M X ,M S ). (37) 

However in minimal SU(5) we find 

e 3 (M x ,Mj,) = 0. (38) 

Thus gauge coupling unification fixes the values of the three parameters, {«g(= 
<7 G /47r), Mq, Mt}- In addition, the ai(Ma) = a2(Mc) condition gives the relation 
( Mg ]1/15 M g M g 

l M T J (M|M S )V3 (^A)i/3(E)- 1 ] 

In the last term we used the approximate relations 

M x ~g G (£), M S ~A(S). (40) 
Hence the natural values for these parameters are given by 

M x ~ M s ~ M G w 3 x 10 16 GeV. (41) 
As a result, the proton lifetime is given by 

^- 10 ^^^n^x!)v r , (42, 



■3 x 10 16 GeV v (3, 



'lattice 
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u <* « u 




Figure 1. X boson exchange diagram giving the dimension 6 four fermion operator 
for proton and neutron decay. 



and the dominant decay mode is 

p -> 7t° + e + . (43) 

Note it is not possible to enhance the decay rate by taking M x <C M G without spoiling 
perturbativity, since this limit requires A ^> 1. On the other hand, Mx S> Mq is 
allowed. 

4-2. Dimension 4 & 5 operators 

The contribution of dimension 4 & 5 operators to nucleon decay in SUSY GUTs was 
noted by [Weinberg (1982); Sakai and Yanagida (1982)]. Dimension 4 operators are 
dangerous. In SUSY GUTs they always appear in the combination 

(U c D c D c ) + (Q L D c ) + (E c L L) (44) 

leading to unacceptable nucleon decay rates. R parity [Farrar and Fayet (1979)] forbids 
all dimension 3 and 4 (and even one dimension 5) baryon and lepton number violating 
operators. It is thus a necessary ingredient of any "natural" SUSY GUT. 

Dimension 5 operators are obtained when integrating out heavy color triplet Higgs 
fields. 

W Z> H u QYJJ + H d (QY d D + LY e E) + (45) 
T{Q l -c qq Q + Uc ue E) + T(Qc q iL + Uc ud D) 

If the color triplet Higgs fields in Eqn. EU{T, T} have an effective mass term My T T 
we obtain the dimension 5 operators 

(1/M^) 



Q-jCqqQ QCqlL + Uc ud D Uc ue E 



(46) 



denoted LLLL and RRRR operators, respectively. Nucleon decay via dimension 5 
operators was considered by [Sakai and Yanagida (1982), Dimopoulos et al (1982), Ellis 
et al (1982)]. The proton decay amplitude is then given generically by the expression 
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c 2 

T(p + u)(x e~FT ( Lo °P Factor) (RG) (K + v\qqql\p) (47) 

M T n 

~ — ^77 (Loop Factor) (RG) ^ lai * tce mp 

M T U J IT 

The last step used a chiral Lagrangian analysis to remove the K + state in favor of the 
vacuum state. Now we only need calculate the matrix element of a three quark operator 
between the proton and vacuum states. This defines the parameter fiiattice using lattice 
QCD calculations. The decay amplitude includes four independent factors: 

(i) Piattice, the three quark matrix element, 

(ii) c 2 , a product of two 3x3 dimensionless coupling constant matrices, 

(iii) a Loop Factor, which depends on the SUSY breaking squark, slepton and gaugino 
masses, and 

(iv) Mj> . , the effective color triplet Higgs mass which is subject to the GUT breaking 
sector of the theory. 

Let us now consider each of these factors in detail. 



4-2.1. f3i a uice The strong interaction matrix element of the relevant three quark 
operators taken between the nucleon and the appropriate pseudo-scalar meson may 
be obtained directly using lattice techniques. However these results have only been 
obtained recently [Aoki et al (JLQCD) (2000), Aoki et al (RBC) (2002)]. Alternatively, 
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chiral Lagrangian techniques [Chadha et al (1983)] are used to replace the pseudo-scalar 
meson by the vacuum. Then the following three quark matrix elements are needed. 

/3C/(k) = e^ 7 <0|(M a rf /3 )M 7 |proton(k)>, (48) 

aU(k) = e Q , / 3 7 <0|(-u* a rf* /3 )-u 7 |proton(k)> (49) 

and C/(k) is the left handed component of the proton's wavefunction. It has been known 
for some time that \/3\ ~ \a\ [Brodsky et al (1984), Gavela et al (1989)] and that \f3\ 
ranges from .003 to .03 GeV 3 [Brodsky et al (1984), Hara et al (1986), Gavela et al 
(1989)]. Until quite recently, lattice calculations did not reduce the uncertainty in \/3\; 
lattice calculations have reported \/3\ as low as .006 GeV 3 [Gavela et al (1989)] and as 
high as .03 GeV 3 [Hara et al (1986)]. Additionally, the phase between a and (3 satisfies 
(3 rs —a [Gavela et al (1989)]. As a consequence, when calculating nucleon decay rates 
most authors have chosen to use a conservative lower bound with |a| ~ \/3\ = 0.003 
GeV 3 and an arbitrary relative phase. 

Recent lattice calculations [Aoki et al (JLQCD) (2000), Aoki et al (RBC) (2002)] 
have obtained significantly improved results. In addition, they have compared the direct 
calculation of the three quark matrix element between the nucleon and pseudo-scalar 
meson to the indirect chiral Lagrangian analysis with the three quark matrix element 
between the nucleon and vacuum. [Aoki et al (JLQCD) (2000)] find 

lattice = (®\qqq\N) = 0.015(1) GeV 3 . (50) 

Also [Aoki et al (RBC) (2002)], in preliminary results reported in conference 
proceedings, obtained 

lattice = 0.007(1) GeV 3 . (51) 

They both find 

^lattice ~ ~ (^lattice- (52) 

Several comments are in order. The previous theoretical range 0.003GeV 3 < Piatuce < 
0.03GeV 3 has been significantly reduced and the relative phase between a and (3 has been 
confirmed. The JLQCD central value is 5 times larger than the previous "conservative 
lower bound." Although the new, preliminary, RBC result is a factor of 2 smaller than 
that of JLQCD. We will have to wait for further results. What about the uncertainties? 
The error bars listed are only statistical. Systematic uncertainties (quenched + chiral 
Lagrangian) are likely to be of order ±50 % (my estimate). This stems from the fact 
that errors due to quenching are characteristically of order 30 %, while the comparison 
of the chiral Lagrangian results to the direct calculation of the decay amplitudes agree 
to within about 20 %, depending on the particular final state meson. 

4-2.2. c 2 - Model Dependence Consider the quark and lepton Yukawa couplings in 
SU(5) - 

A(($)) 10 10 5 H + A'(($)) 10 5 5# (53) 
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or in SO(10) - 

A(($)) 16 16 10 H . (54) 

The Yukawa couplings 

A«$», A'«<t>» (55) 

are effective higher dimensional operators, functions of adjoint (<&) (or higher 
dimensional) representations of SU(5) (or SO(10)). The adjoint representations are 
necessarily there in order to correct the unsuccessful predictions of minimal SU(5) (or 
SO (10)) and to generate a hierarchy of fermion masses. | Once the adjoint (or higher 
dimensional) representations obtain VEVs ((<&)), we find the Higgs Yukawa couplings - 

H u Q Y u U + H d ( Q Y d D + L Y e E) (56) 

and also the effective dimension five operators 

(l/M e T fS ) [ Q X -c qq QQ c ql L + U c ud DU c ue E~Y (57) 

Note, because of the Clebsch relations due to the VEVs of the adjoint 
representations, etc, we have 

7^ C qq 7^ C ue (58) 

and 

Y d ^ Y e ^ c ud ^ c ql . (59) 

Hence, the 3x3 complex matrices entering nucleon decay are not the same 3x3 Yukawa 
matrices entering fermion masses. Is this complication absolutely necessary and how 
large can the difference be? Consider the SU(5) relation - 

\ = K (60) 

[Einhorn and Jones (1982), Inoue et al (1982), Ibanez and Lopez (1984)]. It is known 
to work quite well for small tan (3 ~ 1 — 2 or large tan (3 ~ 50 [Dimopoulos et al (1992), 
Barger et al (1993)]. For a recent discussion see [Barr and Dorsner (2003)]. On the 
other hand, the same relation for the first two families gives 

A s = A M (61) 
Ad = A e 

leading to the unsuccessful relation 

20 ~ ^ = ^ ~ 200. (62) 
m d m e 

This bad relation can be corrected using Higgs multiplets in higher dimensional 
representations [Georgi and Jarlskog (1979), Georgi and Nanopoulos (1979), Harvey et 
al (1980,1982)] or using effective higher dimensional operators [Anderson et al (1994)]. 

J Effective higher dimensional operators may be replaced by Higgs in higher dimensional 
representations, such as 45 of SU(5) or 120 and 126 or SO(10). Using these Higgs representations, 
however, does not by itself address the fermion mass hierarchy. 
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Clearly the corrections to the simple SU(5) relation for Yukawa and c matrices can 
be an order of magnitude. Nevertheless, in predictive SUSY GUTs the c matrices 
are obtained once the fermion masses and mixing angles are fit [Kaplan and Schmaltz 
(1994), Babu and Mohapatra (1995), Lucas and Raby (1996), Frampton and Kong 
(1996), Blazek et al (1997), Barbieri and Hall (1997), Barbieri et al (1997), Allanach 
et al (1997), Berezhiani (1998), Blazek et al (1999,2000), Dermisek and Raby (2000), 
Shafi and Tavartkiladze (2000), Albright and Barr (2000,2001), Altarelli et al (2000), 
Babu et al (2000), Berezhiani and Rossi (2001), Kitano and Mimura (2001), Maekawa 
(2001), King and Ross (2003), Chen and Mahanthappa (2003), Raby (2003), Ross and 
Velasco-Sevilla (2003), Goh et al (2003), Aulakh et al (2003)]. In spite of the above 
cautionary remarks we still find the inexact relations 

c qq Cqi, c ud c ue oc m u m d tan (3. (63) 

In addition, family symmetries can affect the texture of {c qq , c q i, c u d, c ue }, just as it 
will affect the texture of Yukawa matrices. 

In order to make predictions for nucleon decay it is necessary to follow these 
simple steps. Vary the GUT scale parameters, &q, M g , Y u , Y d , Y e and SOFT SUSY 
breaking parameters until one obtains a good fit to the precision electroweak data. 
Whereby we now explicitly include fermion masses and mixing angles in the category of 
precision data. Once these parameters are fit, then in any predictive BGUT the matrices 
c qq , Cqi, c u d, c ue at Mq are also fixed. Now renormalize the dimension 5 baryon and 
lepton number violating operators from Mq — > Mz in the MSSM; evaluate the Loop 
Factor at Mz and renormalize the dimension 6 operators from Mz — > 1 GeV. The latter 
determines the renormalization constant A 3 ~ 1.3 [Dermisek et al (2001)]. [Note, this 
should not be confused with the renormalization factor ~ .22 [Ellis et al (1982)] 
which is used when one does not have a theory for Yukawa matrices. The latter RG 
factor, takes into account the combined renormalization of the dimension 6 operator 
from the weak scale to 1 GeV and also the renormalization of fermion masses from 1 
GeV to the weak scale.] Finally calculate decay amplitudes using a chiral Lagrangian 
approach or direct lattice gauge calculation. 

Before leaving this section we should remark that we have assumed that the 
electroweak Higgs in SO(10) models is contained solely in a 10. If however the 
electroweak Higgs is a mixture of weak doublets from 16 h, 10 h and, in addition, we 
include the higher dimensional operator jj(16 16 16# 16#), useful for neutrino masses, 
then there are additional contributions to the dimension 5 operators considered in (Eqn. 
l4*oT) [Babu et al (2000)]. However these additional terms are not required for neutrino 
masses [Blazek et al (1999,2000)]. 

4 . 3. Loop factor 

The dimension 5 operators are a product of two fermion and two scalar fields. The scalar 
squarks and/or sleptons must be integrated out of the theory below the Bbreaking scale. 
There is no consensus on the best choice for an appropriate Bbreaking scale. Moreover, 
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Figure 3. Higgsino and third generation stop - stau loop giving the dominant loop 
contribution to dimension 5 nucleon decay. 



in many cases there is a hierarchy of fiparticle masses. Hence we take the simplest 
assumption, integrating out all fiparticles at M z . When integrating out the Bparticles 
in loops, the effective dimension 5 operators are converted to effective dimension 6 
operators. This results in a loop factor which depends on the sparticle masses and 
mixing angles. 

Consider the contribution to the loop factor for the process p — > K + + v T in Fig. El 
This graph is due to the RRRR operators and gives the dominant contribution at large 
tan (5 and a significant contribution for all values of tan (3 [Lucas and Raby (1997), Babu 
and Strassler (1998), Goto and Nihei (1999), Murayama and Pierce (2002)]. Although 
the loop factor is a complicated function of the sparticle masses and mixing angles, 
it nevertheless has the following simple dependence on the overall gaugino and scalar 
masses given by 



W a+ (64) 

167T 2 mf 6 
Thus in order to minimize this factor one needs 

(i, Mi /2 SMALL (65) 

and 



mie Large. 



(66) 
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The largest uncertainty in the nucleon decay rate is due to the color triplet Higgs mass 
parameter M T . As M T increases, the nucleon lifetime increases. Thus it is useful 
to obtain an upper bound on the value of My . This constraint comes from imposing 
perturbative gauge coupling unification [Lucas and Raby (1997), Goto and Nihei (1999), 
Babu et al (2000), Altarelli et al (2000), Dermfsek et al (2001), Murayama and Pierce 
(2002)]. Recall, in order to fit the low energy data a GUT scale threshold correction 

e 3 ss MMg) " &G) ~ -3% to - 4% (67) 

is needed. €3 is a logarithmic function of particle masses of order Mq, with contributions 
from the electroweak Higgs and GUT breaking sectors of the theory. 

e 3 = e f ggs + e£ UT breaking + • • • , (68) 

^If'-C'' (69) 

In TableEJwe have analyzed three different GUT theories - the minimal SU(5) model, an 
SU(5) model with natural doublet-triplet splitting and minimal SO(10) (which also has 
natural doublet-triplet splitting). We have assumed that the low energy data, including 
weak scale threshold corrections, requires e 3 = —4%. We have then calculated the 
contribution to €3 from the GUT breaking sector of the theory in each case. 

Minimal SU(5) is defined by its minimal GUT breaking sector with one SU(5) 
adjoint S. The one loop contribution from this sector to £3 vanishes. Hence, the 
—4 % must come from the Higgs sector alone, requiring the color triplet Higgs mass 
~ 2 x 10 14 GeV. Note since the Higgs sector is also minimal, with the doublet 
masses fine-tuned to vanish, we have M T = My. By varying the Bspectrum at the 
weak scale, we may be able to increase €3 to —3 % or even —2 %, but this cannot save 
minimal SU(5) from disaster due to rapid proton decay from Higgsino exchange. 

In the other theories, Higgs doublet - triplet splitting is obtained without fine- 
tuning. This has two significant consequences. First, the GUT breaking sectors are 
more complicated, leading in these theories to large negative contributions to €3. The 
maximum value | e 3 1 ~ 10% in minimal SO(10) is fixed by perturbativity bounds 
[Dermisek et al (2001)]. Secondly, the effective color triplet Higgs mass M^ ff does 
not correspond to the mass of any particle in the theory. In fact, in both cases with 
"natural" doublet-triplet splitting, the color triplet Higgs mass is of order M G even 
though M T 3> Mq. The values for MJ in Table El are fixed by the value of e|j? lggh 
needed to obtain e 3 = —4%. 

Before discussing the bounds on the proton lifetime due to the exchange of color 
triplet Higgsinos, let us elaborate on the meaning of My . Consider a simple case with 
two pairs of SU(5) Higgs multiplets, { 5 l H , 5 l H } with i = 1,2. In addition, also assume 
that only {5^, 5^ } couples to quarks and leptons. Then My is defined by the 
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Table 3. GUT threshold corrections in three different theories. The first column 
is the minimal SU(5) fitheory [Dimopoulos and Georgi (1981) Sakai (1982), Goto and 
Nihei (1999), Murayama and Pierce (2002)], the second column is SU(5) with "natural" 
Higgs doublet-triplet splitting [Masiero et al (1982), Grinstein (1982), Altarelli et al 
(2000)], and the third column is minimal SO(10) fimodel [Blazek et al (1999.2000), 
Dermisek et al (2001), Blazek et al (2002a,b)]. 





Minimal SU(5) 


SU(5) "Natural" D/T 


Minimal SO(10) 


GUTbrcaking 
£ 3 





- 7.7 % 


- 10 % 


Higgs 
e 3 


- 4 % 


+ 3.7 % 


+ 6 % 


M* ff [GeV] 


2 x 10 14 


3 x 10 18 


6 x 10 19 



expression 

^ = <«?% (™) 

where Mt is the color triplet Higgs mass matrix. In the cases with "natural" doublet - 
triplet splitting, we have 

Ur-(.°. M °) PI) 



X 



with 



X 



M T ff ~M 2 G (72) 

Thus for X <C Mq we have M T ^ 3> Mq and no particle has mass greater than Mq 
[Babu and Barr (1993)]. The large Higgs contribution to the GUT threshold correction 
is in fact due to an extra pair of light Higgs doublets with mass of order X. 

Due to the light color triplet Higgsino, it has been shown that minimal SUSY 
SU(5) is ruled out by the combination of proton decay constrained by gauge coupling 
unification [Goto and Nihei (1999), Murayama and Pierce (2002)] !! In Figs. 0] and 
13 we reprint the figures from the paper by [Goto and Nihei (1999)]. In Fig. Hj the 
decay rate for p — > K + j?i for any one of the three neutrinos (i = e,/i and r) is 
plotted for fixed soft Bbreaking parameters as a function of the relative phase 023 
between two LLLL contributions to the decay amplitude. The phase 0i3 is the relative 
phase between one of the LLLL contributions and the RRRR contribution. The latter 
contributes predominantly to the v T final state, since it is proportional to the up quark 
and charged lepton Yukawa couplings. As noted by [Goto and Nihei (1999)], the partial 
cancellation between LLLL contributions to the decay rate is completely filled by the 
RRRR contribution. It is this result which provides the stringent limit on minimal 
fiSU(5). As one sees from Fig. HI for the color triplet Higgs mass Mt = 2 x 10 16 GeV 
(= Mc in the notation of [Goto and Nihei (1999)]), the universal scalar mass uiq = 
1 TeV and tan/3 = 2.5, there is no value of the phase 023 which is consistent with 
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Figure 4. [Fig. 2: Goto and Nihei (1999)] Decay rates T(p -> #+£>,) {i = e,n and r) 
as functions of the phase </>23 for tan/3 = 2.5. The soft Bbreaking parameters are 
fixed at mo = 1 TcV, M- g = 125 GeV and A = with /u > 0. M T and M s are 
taken as Mt = Ms = 2 x 10 16 GeV. The horizontal lower line corresponds to the 
Super-Kamiokande limit r(p — > > 5.5 x 10 32 years, and the horizontal upper 

line corresponds to the Kamiokande limit r(p — > K + v) > 1.0 x 10 32 years. The new 
Super-Kamiokande bound is r(p — > i^T+P) > 1.9 x 10 33 years. 



Super-Kamiokande bounds. Note, the proton decay rate scales as tan/3 2 ; hence the 
disagreement with data only gets worse as tan /3 increases. In Fig. EJ the contour of 
constant proton lifetime is plotted in the Mt(= Mc) - m(ui) plane, where m(ui) is the 
mass of the left-handed up squark for tan/? = 2.5. Again, there is no value of m(ui) < 3 
TeV for which the color triplet Higgs mass is consistent with gauge coupling unification. 
In [Goto and Nihei (1999)] the up squark mass was increased by increasing m . Hence 
all squarks and slepton masses increased. 
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Figure 5. [Fig. 4: Goto and Nihei (1999)] Lower bound on the colored Higgs mass 
Mt as a function of the left-handed scalar up quark mass m^ L ■ The soft fibrcaking 
parameters m , Mg, A are scanned within the range < m < 3 TeV, < Mg < 1 
TeV and — 5 < Ao < 5 with tan/3 fixed at 2.5. Both signs of fi are considered. 
The whole parameter region of the two phases <f>ia and <p23 is examined. The solid 
curve represents the bound derived from the Super-Kamiokande limit r(p — > K + v) > 
5.5 x 10 32 years and the dashed curve represents the corresponding result without 
the RRRR effect. The left-hand side of the vertical dotted line is excluded by other 
experimental constraints. The dash-dotted curve represents the bound derived from 
the Kamiokande limit on the neutron partial lifetime r(n — > K°D) > 0.86 x 10 32 years. 



One may ask whether one can suppress the proton decay rate by increasing the 
mass of the squarks and sleptons of the first and second generation, while keeping the 
third generation squarks and sleptons light (in order to preserve "naturalness"). This 
is the question addressed by [Murayama and Pierce (2002)]. They took the first and 
second generation scalar masses of order 10 TeV, with the third generation scalar masses 
less than 1 TeV. They showed that since the RRRR contribution does not decouple in 
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this limit, and moreover since any possible cancellation between the LLLL and RRRR 
diagrams vanishes in this limit, one finds that minimal 6SU(5) cannot be saved by 
decoupling the first two generations of squarks and sleptons. 

Thus minimal BSU(5) is dead. Is this something we should be concerned about. In 
my opinion, the answer is no, although others may disagree [Bajc et al (2002)]. Minimal 
BSU(5) has two a priori unsatisfactory features: 

• It requires fine-tuning for Higgs doublet-triplet splitting, and 

• renormalizable Yukawa couplings due to 5#, 5# alone are not consistent with 
fermion masses and mixing angles. 

Thus it was clear from the beginning that two crucial ingredients of a realistic theory 
were missing. The theories which work much better have "natural" doublet-triplet 
splitting and fit fermion masses and mixing angles. 



4-5. Summary of Nucleon Decay in J^D 

Minimal BSU(5) is excluded by the concordance of experimental bounds on proton decay 
and gauge coupling unification. We discussed the different factors entering the proton 
decay amplitude due to dimension 5 operators. 



T{p ^K + + v)~ -^(Loop Factor) =^ m p . (73) 

We find 

• c 2 : model dependent but constrained by fermion masses and mixing angles; 

• ^lattice' JLQCD central value is 5 times larger than the previous "conservative 
lower bound." However one still needs to reduce the systematic uncertainties of 
quenching and chiral Lagrangian analyses. Moreover, the new RBC result is a 
factor of 2 smaller than JLQCD; 



• Loop Factor: oc x * 6 ^£ - m ^ — — . It is minimized by taking gauginos light and the 1st 
and 2nd generation squarks and sleptons heavy (> TeV). However, "naturalness" 
requires that the stop, sbottom and stau masses remain less than of order 1 TeV; 

• . constrained by gauge coupling unification and GUT breaking sectors. 

The bottom line we find for dimension 6 operators [Lucas and Raby (1997), 
Murayama and Pierce (2002)] 

+x r n n36/ M x x4 ,0.015 GeV 3 . 2 
r (P ^' + ( >5x 10- ( 3xl016GeV )' l-j— f <") 

Note, it has been recently shown [Klebanov and Witten (2003)] that string theory can 
possibly provide a small enhancement of the dimension 6 operators. Unfortunately the 
enhancement is very small. Thus it is very unlikely that these dimension 6 decay modes 
p — > 7T° + e + will be observed. 
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On the other hand for dimension 5 operators in realistic SUSY GUTs we obtain 
rough upper bounds on the proton lifetime coming from gauge coupling unification and 
perturbativity [Babu et al (2000), Altarelli et al (2000), Dermisek et al (2001)] 

T (p _> K + + v) < (\ - 3) x 10 34 ( ° M5 GeV V years. (75) 

<J Plattice 

Note in general 

T(n-*K° + v) <t(p^K + + v). (76) 

Moreover other decay modes may be significant, but they are very model dependent, for 
example [Carone et al (1996), Babu et al (2000)] 

p-^7r° + e + , K° + /i + . (77) 



4-6. Proton decay in more than four dimensions 

We should mention that there has been a recent flurry of activity on SUSY GUTs in extra 
dimensions beginning with the work of [Kawamura (2001a,b)]. However the study of 
extra dimensions on orbifolds goes back to the original work of [Dixon et al (1985,1986)] 
in string theory. Although this interesting topic would require another review, let me 
just mention some pertinent features here. In these scenarios, grand unification is only 
a symmetry in extra dimensions which are then compactified at scales of order 1/Mq. 
The effective four dimensional theory, obtained by orbifolding the extra dimensions, has 
only the standard model gauge symmetry or at most a Pati-Salam symmetry which 
is then broken by the standard Higgs mechanism. In these theories, it is possible to 
completely eliminate the contribution of dimension 5 operators to nucleon decay. This 
may be a consequence of global symmetries as shown by [Witten (2002), Dine et al 
(2002)] or a continuous U(1)_r symmetry (with R parity a discrete subgroup) [Hall and 
Nomura (2002)]. [Note, it is also possible to eliminate the contribution of dimension 
5 operators in 4 dimensional theories with extra symmetries [Babu and Barr (2002)], 
but these 4 dimensional theories are quite convoluted. Thus it is difficult to imagine 
that nature takes this route. On the other hand, in one or more small extra dimensions 
the elimination of dimension 5 operators is very natural.] Thus at first glance, nucleon 
decay in these theories may be extremely difficult to see. However this is not necessarily 
the case. Once again we must consider the consequences of grand unification in extra 
dimensions and gauge coupling unification. 

Extra dimensional theories are non-renormalizable and therefore require an explicit 
cutoff scale M*, assumed to be larger than the compactification scale M c . The Kaluza- 
Klein excitations above M c contribute to threshold corrections to gauge coupling 
unification evaluated at the compactification scale. The one loop renormalization of 
gauge couplings is given by 

2tt 2tt ,M„. 
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where A, are the threshold corrections due to all the KK modes from M c to M* and can 
be expressed as A< = bf f log(f|) [Hall and Nomura (2001,2002a), Nomura et al (2001), 
Contino et al (2002), Nomura C (2002)]. In a 5D SO(10) model, broken to Pati-Salam 
by orbifolding and to the MSSM via Higgs VEVs on the brane, it was shown [Kim and 
Raby (2003)] that the KK threshold corrections take a particularly simple form 

2 " 2 " +6f 5SM log(^)+A l (79) 



Oii(fi) a(M*) ' n 

2 71 " , iMSSMi , ^uSM/ir\i /M* 



+ bf log(— ) + =-bt M (V) log(^). (80) 



with 



Agauge =3 & rG0log(^) (81) 



^Higgs — 0, (82) 

and M c = (nR/2)^ 1 . Here h MSSM includes the gauge sector and the Higgs sector 
together and b (V) includes the gauge sector only. The running equation is very 
simple and permits us to directly compare with well known 4D SUSY GUTs. In the 
minimal 4D SU(5) Bmodel, the running equation is given by 

— pr = + b? SSM log(^) + 6f M ( V) log(^). (83) 

OLi(n) a{M GUT ) /i M T 

where Mr is the color triplet Higgs mass. As discussed earlier, we can achieve unification 
by adjusting the color triplet Higgs mass M T = 2 x 10 14 GeV (see Table |3J). 
Comparing Eqn. (|50|) and (|S3^1. we observe that if 

M C = M T , ^— = ( )2 

M c Mt 

the same unification is achieved here. Therefore we get M c = 2 x 10 14 GeV and 
M* = 3.7 x 10 17 GeV by using Mqut = 3 x 10 16 GeV. A few remarks are in order. 
There is no problem with proton decay due to dimension 5 operators, even though 
the color triplet Higgs mass is of order 10 14 GeV, since these operators are excluded 
by R symmetry. In a 5D model, the 4D GUT scale has no fundamental significance. 
The couplings unify at the cutoff scale and there is no scale, above which we have a 
perturbative SO(10) GUT. 

In Fig. El we show the running in the difference of couplings for two independent 
cases. 

(i) We show the couplings for four dimensional gauge theories with GUT scale 
thresholds, in which the GUT scale is defined as the point where ol\ and «2 meet 
and £3 is the relative shift in due to threshold corrections, and 

(ii) for a five dimensional SO (10) model where the three couplings meet at the 
cutoff scale and the threshold corrections due to the KK tower is defined at the 
compactification scale. 
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Figure 6. Differential running of do — 2tt(— — ) and 63 = 2n(— —). 

In both cases, the running of the gauge couplings below the compactification scale must 
be the same. Thus we can use the low energy fits from 4D theories to constrain a 5D 
theory. 

Note, since the the KK modes of the baryon number violating {X, X} gauge bosons 
have mass starting at the compactification scale M c « 10 14 GeV we must worry whether 
proton decay due to dimension 6 operators is safe. It has been shown [Altarelli and 
Feruglio (2001), Hall and Nomura (2002b), Hebecker and March-Russell (2002)] that 
this depends on where in the extra dimensions the quarks and leptons reside. If they 
are on symmetric orbifold fixed points, i.e. symmetric under the GUT symmetry, then 
this leads to the standard dimension 6 proton decay operators which is ruled out for 
the first and second families of quarks and leptons. Hence the first two families must 
be either in the bulk or on broken symmetry fixed points. If they are in the bulk, then 
the {X, X} mediated processes take massless modes to KK excitations which is not a 
problem. Otherwise, if the first two families are on the broken symmetry fixed point, 
the wave functions for the {X, X} bosons vanish there. However, certain effective 
higher dimensional operators on the broken symmetry orbifold fixed points can allow 
the {X, X} bosons to couple to the first two families. These operators are allowed 
by symmetries and they naturally lead to proton lifetimes for p — > e + tt° of order 
10 34±2 years. The large uncertainty is due to the order of magnitude uncertainty in the 
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coefficient of these new effective operators. 

Before closing this section, we should make a few comments on theories with large 
extra dimensions of order 1/TeV or as large as 1 mm. In some of these theories, only 
gravity lives in higher dimensions while the ordinary matter and gauge interactions 
typically reside on a three dimensional brane [Arkani-Hamed et al (1998)]. While if 
the extra dimensions are no larger than 1/TeV, all matter may live in the bulk. Such 
theories replace SUSY with new and fundamental non-perturbative physics at the 1-10 
TeV scale. These theories must address the question of why dimension 6 proton decay 
operators, suppressed only by 1/TeV 2 , are not generated. There are several ideas in the 
literature with suggested resolutions to this problem. They include: 

• conserved baryon number on the brane with anomaly cancelling Chern-Simons 
terms in the bulk, 

• or displacing quarks from leptons on a "fat" brane with a gaussian suppression of 
the overlap of the quark/lepton wave functions. 

A novel solution to the problem of proton decay is found in 6 space-time dimensional 
theories with all matter spanning the two extra dimensions. It has been shown that 
in such theories [Appelquist et al (2001)] a Z(8) remnant of the 6 dimensional Lorentz 
group is sufficient to suppress proton decay to acceptable levels. These theories predict 
very high dimensional proton decay operators with multi-particle final states. 

5. Fermion masses and mixing 

Low energy Bprovides a natural framework for solving the gauge hierarchy problem, 
while SUSY GUTs make the successful prediction for gauge coupling unification and 
the, still un-verified, prediction for proton decay. But these successes affect only a small 
subset of the unexplained arbitrary parameters in the standard model having to do with 
the Z and Higgs masses (i.e. the weak scale) and the three gauge coupling constants. On 
the other hand, the sector of the standard model with the largest number of arbitrary 
parameters has to do with fermion masses and mixing angles. Grand unification also 
provides a natural framework for discussing the problem of fermion masses, since it 
naturally arranges quarks and leptons into a few irreducible multiplets, thus explaining 
their peculiar pattern of gauge charges, i.e. charge quantization and the family structure. 
Moreover, it has been realized for some time that the masses and mixing angles of quarks 
and leptons are ordered with respect to their generation (or family) number. The first 
generation of quarks and leptons, {u, d, e, ^ e }, are the lightest; the second generation, 
{c, s, [A, u^}, are all heavier than the first, and the third generation, {£, b, r, i/ r },arethe 
heaviest (see Table EJ). In addition, the first two generations have a weak mixing angle 
given by the Cabibbo angle, 9c- If we define A = sin^c ~ -22, then the mixing of the 
second and third generation is of order A 2 and the first and third is the weakest mixing 
of order A 3 . This pattern is very elegantly captured in the Wolfenstein representation 
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Table 4. Quark and lepton masses in units of MeV/c 2 . 







e 


u 


d 


1st 


< 10- 7 


1/2 


2 


5 


2nd 


< 10" 7 


105.6 


1,300 


120 


3rd 


< 10~ 7 


1,777 


174,000 


4,500 



of the CKM matrix given by 

/ V ud 

VcKM = Vcd 

V v td 

I 



v us 
v cs 

Vt 



ts 



v cb 
v tb ) 

\* 

-X 

irj) A A 3 
~ 1/2. 



M) 



A 

-A A 2 



(p - irf) A A 3 \ 
A A 2 

1 J 



with A ~ 1/5, A ~ 1 and |p — zry| 

Although the fundamental explanation for three families is still wanting, it is natural 
to assume that the families transform under some family symmetry. Such a possibility is 
consistent with weakly coupled heterotic string theory where, for example, E(8) x E(8) 
and SO (32) are both large enough to contain a GUT group x a family symmetry. Such 
a family symmetry has many potential virtues. 

• A spontaneously broken family symmetry can explain the hierarchy of fermion 
masses and mixing angles. 

• In a Btheory, the family symmetry acts on both fermions and sfermions, thus 
aligning the quark and squark, and the lepton and slepton mass matrices. This 
suppresses flavor violating processes. 

• The combination of a family and a GUT symmetry can reduce the number of 
fundamental parameters in the theory, hence allowing for a predictive theory. 

In the following sections, we consider several important issues. In section 15.11 we 
discuss the simplest case of the third generation only. Here we discuss the status of 
SU(5) (A 6 = A T ) and SO(10) (X t = X b = X T = X Vt ) Yukawa unification. In section 
15.21 we consider several different analyses in the literature for three family models with 
either U(l) or non-abelian family symmetry. In section l5~3*l we study the relation between 
charged fermion and neutrino masses in SUSY GUTs and consider some examples giving 
bi-large neutrino mixing consistent with the data. Finally, in section 15^1 we discuss some 
experimental consequences of fitheories of fermion masses. In particular, we consider 



si, (9 



2) M , B s 



ji + fi ,[/,—> e^y and the electric dipole moments d n , d e . 
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5.1. Yukawa unification 

Let us first discuss the most stringent case of SO (10) Yukawa unification. It has been 
shown [Raby (2001), Derrmsek (2001), Baer and Ferrandis (2001), Blazek et al (2002a,b), 
Auto et al (2003), Tobe and Wells (2003)] that SO(10) boundary conditions at the GUT 
scale, for soft SUSY breaking parameters as well as for the Yukawa couplings of the third 
generation, are consistent with the low energy data, including M t , m b (mb), M T , ONLY 
in a narrow region of SUSY breaking parameter space. Moreover, this region is also 
preferred by constraints from CP and flavor violation, as well as by the non-observation 
of proton decay. Finally we discuss the consequences for the Higgs and SUSY spectrum. 

Recall, in SO(10) we have the compelling unification of all quarks and leptons 
of one family into one irreducible representation such that 10 + 5 + v steri i e C 16 
and the two Higgs doublets are also unified with 5h, 5h C IOh- Hence, minimal 
SO(10) also predicts Yukawa unification for the third family of quarks and leptons 
with Afe = A* = A T = \ Ut = A at the GUT scale [Banks (1988), Olechowski and 
Pokorski (1988), Pokorski (1990), Shafi and Ananthanarayan (1991);Ananthanarayan et 
al (1991,1993,1994), Anderson et al (1993,1994)]. 

Ignoring threshold corrections, one can use the low energy value for mj/m T to fix 
the universal Yukawa coupling A. RG running from Mq to Mz then gives X T (Mz)- Then 
given m T = \ T -^cosf3 we obtain tan/5 ~ 50. Finally, a prediction for the top quark 
mass is given with m t = X t ^sin/3 ~ 170 ± 20 GeV (see [Anderson et al (1993)]). 

Note, in this case there are insignificant GUT threshold corrections from gauge and 
Higgs loops. Nevertheless, the previous discussion is essentially a straw man, since there 
are huge threshold corrections at the weak scale [Hall et al (1994), Hempfling (1994), 
Carena et al (1994), Blazek et al (1995)]. The dominant contributions are from gluino 
and chargino loops plus an overall logarithmic contribution due to finite wave function 
renormalization given by bm^jm^ = Am\ + Am% + Am\ og + • • • (see Fig. EJ). These 
contributions are approximately of the form 

Am»^^ tan(3, (85) 
37r m{ 

b 

Am? « ^ tanp and (86) 

V6i\ l m~ 

Am^glog(^-)~6% (87) 
z 

with Am 9 b ~ -AmJ > for fi > [with our conventions] . These corrections can easily 
be of order ~ 50 %. However good fits require —4% < 5m b /m b < —2%. 

Note, the data favors \x > 0. First consider the process b — > 57. The chargino 
loop contribution typically dominates and has opposite sign to the standard model 
and charged Higgs contributions for /i > 0, thus reducing the branching ratio. This 
is desirable, since the standard model contribution is a little too large. Hence \x < 
is problematic when trying to fit the data. Secondly, the recent measurement of the 
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Figure 7. The one loop gluino (left) and chargino (right) corrections to the bottom 
quark mass proportional to a s (left) and X t (right) and to tan (3. 



anomalous magnetic moment of the muon suggests a contribution due to NEW physics 
given by a^ EW = 22.1(11.3) x 1(T 10 or 7.4(10.5) x 10~ 10 [Muon g - 2 Collaboration 
(2002), Davier et al (2003)] depending on whether one uses e + e _ or r hadronic decay 
data to evaluate the leading order hadronic contributions. For other recent theoretical 
analyses and references to previous work see [Hagiwara et al (2003), Melnikov and 
Vainshtein (2003)]. However in SUSY the sign of a^ EW is correlated with the sign of \x 
[Chattopadhyay and Nath (1996)]. Once again the data favors \x > 0. 

Before discussing the analysis of Yukawa unification, specifically that of [Blazek 
et al (2002a,b)], we need to consider one important point. 50(10) Yukawa unification 
with the minimal Higgs sector necessarily predicts large tan/5 ~ 50. In addition, it 
is much easier to obtain EWSB with large tan/3 when the Higgs up/down masses are 
split {tti 2 h u < m // d ) [Olechowski and Pokorski (1995), Matalliotakis and Nilles (1995), 
Polonsky and Pomarol (1995), Murayama, Olechowski and Pokorski (1996), Rattazzi 
and Sarid (1996)]. In the following analysis we consider two particular Higgs splitting 
schemes we refer to as Just So and D term splitting. § In the first case the third generation 
squark and slepton soft masses are given by the universal mass parameter m 16 , and only 
Higgs masses are split: Tn 2 Hu H ^ = mf (1 qp Am#). In the second case we assume D 
term splitting, i.e. that the D term for U(l)x is non-zero, where U(l)x is obtained 
in the decomposition of 50(10) — > SU(b) x U(l)x- In this second case, we have 
m h u ,H d ) = m lo T 2D X , m(Q,u,e) = m\ & + D x , m 2 (lL) = m 2 w - 3D X . The Just So 
case does not at first sight appear to be very well motivated. However we now argue 
that it is quite natural [Blazek et al (2002a,b)]. In 5O(10), neutrinos necessarily have 
a Yukawa term coupling active neutrinos to the "sterile" neutrinos present in the 16. 
In fact for u T we have L Ut v t L H u with ~L Ut = L t = L& = L T = L. In order to 
obtain a tau neutrino with mass vh Vt ~ 0.05 eV (consistent with atmospheric neutrino 
oscillations), the "sterile" v T must obtain a Majorana mass M Pr > 10 13 GeV. Moreover, 
since neutrinos couple to H u (and not to Hd) with a fairly large Yukawa coupling (of order 

§ Just So Higgs splitting has also been referred to as non universal Higgs mass splitting or NUHM 
[Berezinsky et al (1996), Blazek et al (1997a,b), Nath and Arnowitt (1997)]. 
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Figure 8. x 2 contours for m 1& — 1.5 TeV (Left) and mi 6 = 2 TeV (Right). The 
shaded region is excluded by the chargino mass limit > 103 GeV. 

0.7), they naturally distinguish the two Higgs multiplets. With L = 0.7 and M Vt = 10 14 
GeV, we obtain a significant GUT scale threshold correction with Am 2 H ~ 7%, about 
1/2 the value needed to fit the data. At the same time, we obtain a small threshold 
correction to Yukawa unification m 1.75%. 

5.1.1. x 2 Analysis [Blazek et al (2002a, b)] Our analysis is a top-down approach 
with 11 input parameters, defined at Mq, varied to minimize a x 2 function composed 
of 9 low energy observables. The 11 input parameters are: Mg, ocg(Mg), 63; the 
Yukawa coupling L, and the 7 soft SUSY breaking parameters /1, M1/2, A , tan/3, 
m? 6 , ml , Am 2 H (Dx) for Just So (D term) case. We use two (one)loop renormalization 
group [RG] running for dimensionless (dimensionful) parameters from Mq to Mz and 
complete one loop threshold corrections at Mz [Pierce et al (1997)]. We require 
electroweak symmetry breaking using an improved Higgs potential, including m\ and 
mf corrections in an effective 2 Higgs doublet model below M stop [Haber and Hempfling 
(1993), Carena et al (1995,1996)]. Note, in the figures we have chosen to keep three 
input parameters /1, M1/2, mie fixed, minimizing x 2 with respect to the remaining 8 
parameters only. The \ 2 function includes the 9 observables; 6 precision electroweak 
data «em, Gfi, oc s (Mz) = 0.118 (0.002), M z , M w , Pnew an d the 3 fermion masses 
M top = 174.3 (5.1), m b (m b ) = 4.20 (0.20), M T . Fig. H (Left) shows the constant \ 2 
contours for mie = 1.5 TeV in the case of Just So squark and slepton masses. We find 
acceptable fits (x 2 < 3) for A ~ —1.9 mie, m w ~ 1.4 mie and mie > 1.2 TeV. The best 
fits are for > 2 TeV with x 2 < 1. Fig. 1 (Right) shows the constant \ 2 contours for 
mi6 = 2 TeV. Fig. Ogives the constant m^mf,) and bm^jm^ contours for mig = 2 TeV. 
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Figure 9. Contours of constant mb(mb)[GeV] (Left) and Amj in % (Right) for 
mie = 2 TeV. 

We see that the best fits, near the central value, are found with —4% < 5m b /m b < —2%. 
The chargino contribution (Eqn. IHUj) is typically opposite in sign to the gluino (Eqn. 

since A t runs to an infrared fixed point oc — Mi/ 2 (see for example, [Carena e£ a/ 
(1994)]). Hence in order to cancel the positive contribution of both the log (Eqn. IHTj) 
and gluino contributions, a large negative chargino contribution is needed. This can 
be accomplished for — A t > rrig and <C mg . The first condition can be satisfied 
for Aq large and negative, which helps pull A t away from its infrared fixed point. The 
second condition is also aided by large A t . However in order to obtain a large enough 
splitting between and m~ bi , large values of are needed. Note, that for Just So 
scalar masses, the lightest stop is typically lighter than the sbottom. We typically find 
mg ~ 3 mf x . On the other hand, D term splitting with Dx > gives mg < m^. As 
a result in the case of Just So boundary conditions excellent fits are obtained for top, 
bottom and tau masses; while for D term splitting the best fits give m^m^) > 4.59 
GeV.|| 

|| Note, [Auto et al (2003), Tobe and Wells (2003)] use a bottom-up approach in their analysis. The 
results of [Auto et al (2003)] are in significant agreement with [Blazek et al (2002a, b)], except for the 
fact that they only find Yukawa unification for larger values of m\Q of order 8 TeV and higher. The 
likely reason for this discrepancy has been explained by [Tobe and Wells (2003)]. They show that 
[and I quote them] "Yukawa couplings at the GUT scale arc very sensitive to the low-energy SUSY 
corrections. An 0(1%) correction at low energies can generate close to a O(10%) correction at the 
GUT scale. This extreme IR sensitivity is one source of the variance in conclusions in the literature. 
For example, course-grained scatter plot methods, which are so useful in other circumstances, lose some 
of their utility when IR sensitivity is so high. Furthermore, analyses that use only central values of 
measured fermion masses do not give a full picture of what range of supersymmetry parameter space 
enables third family Yukawa unification, since small deviations in low-scale parameters can mean so 
much to the high-scale theory viability." It should also be noted that [Tobe and Wells (2003)] suggest 
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The bottom line is that Yukawa unification is only possible in a narrow region of 
SUSY parameter space with 

A ~ —1.9 mie, mio ~ 1.4 mie, (88) 
(//, M V2 ) ~ 100 - 500 GeV and m 16 > 1.2 TeV. (89) 



It would be nice to have some a priori reason for the fundamental SUSY breaking 
mechanism to give these soft SUSY breaking parameters. However, without such an 
a priori explanation, it is all the more interesting and encouraging to recognize two 
additional reasons for wanting to be in this narrow region of parameter space. 

(i) One mechanism for suppressing large flavor violating processes in SUSY theories is 
to demand heavy first and second generation squarks and sleptons (with mass 3> 
TeV) and the third generation scalars lighter than a TeV. Since the third generation 
scalars couple most strongly to the Higgs, this limit can still leave a "naturally" 
light Higgs [Dimopoulos and Giudice (1995)]. It was shown that this inverted 
scalar mass hierarchy can be obtained "naturally," i.e. purely as a consequence of 
renormalization group running from Mq to Mz, with suitably chosen soft SUSY 
breaking boundary conditions at M G [Bagger et al (1999,2000)]. All that is needed 
is SO (10) boundary conditions for the Higgs mass (i.e. mio), squark and slepton 
masses (i.e. m^) an d a universal scalar coupling A . In addition, they must be in 
the ratio [Bagger et al (2000)] 

Al = 2m 2 w = A mj e , with m 16 > TeV. (90) 

(ii) In order to suppress the rate for proton decay due to dimension 5 operators one 
must also demand [Dermisek et al (2001)] 

(fi, Mi/ 2 ) <C mi 6 , with mi 6 > few TeV. (91) 

5.1.2. Consequences for Higgs and SUSY Searches In Fig. El we show the constant 
light Higgs mass contours for m 16 = 1.5 and 2 TeV (solid lines) with the constant y 2 
contours overlayed (dotted lines). Yukawa unification for \ 2 < 1 clearly prefers a light 
Higgs with mass in a narrow range, 112 - 118 GeV. 

In this region the CP odd A, the heavy CP even Higgs H and the charged Higgs 
bosons are also quite light. In addition we find the mass of ti ~ (150 — 250) GeV, 
bi ~ (450 - 650) GeV, h ~ (200 - 500) GeV, g ~ (600 - 1200) GeV, x + ~ (100 - 250) 
GeV, and x° ~ (80 — 170) GeV. All first and second generation squarks and sleptons 
have mass of order m 16 . The light stop and chargino may be visible at the Tevatron. 
With this spectrum we expect ti — » x + b with x + ~ > Xi I v 1° t" e dominant. Lastly Xx 
is the LSP and possibly a good dark matter candidate (see for example, [Roszkowski et 
al (2001)] and Fig. EJ. 

a different soft fibreaking solution consistent with Yukawa unification. In particular, they suggest an 
extension of AMSB with the addition of a large universal scalar mass mo > 2 TeV. 



Desperately Seeking Supersymmetry [SUSY] 



44 




Figure 10. Contours of constant rrih [GeV] (solid lines) with x 2 contours from Fig. 
1 (dotted lines) for m 16 = 1.5 TeV (Left) and m 16 = 2 TeV (Right). 



Our analysis thus far has only included third generation Yukawa couplings; hence 
no flavor mixing. If we now include the second family and 2-3 family mixing, consistent 
with V c b, we obtain new and significant constraints on mj and m^. The stop mass is 
constrained by B(b — > 57) to satisfy m~ IIN > 450 GeV (unfortunately increasing the 
bottom quark mass). In addition, as shown by [Choudhury and Gaur (1999), Babu 
and Kolda (2000), Dedes et al (2001), Isidori and Retico (2001)] the one loop SUSY 
corrections to CKM mixing angles (see Blazek et al (1995)) result in flavor violating 
neutral Higgs couplings. As a consequence the CDF bound on the process B s — > ji + 
places a lower bound on uia > 200 GeV [Choudhury and Gaur (1999), Babu and Kolda 
(2000), Dedes et al (2001), Isidori and Retico (2001)]. x 1 , on the other hand, increases 
as m^o increases. However the increase in y 2 is less than 60% for < 400 GeV. Note, 
the H , H° masses increase linearly with itla- 

5.1.3. SU(5) Yukawa unification Now consider Yukawa unification in SU(5). In this 
case we only have the GUT relation A& = A r . The RG running of the ratio \b/\ to 
low energies then increases (decreases) due to QCD (Yukawa) interactions. In addition, 
neglecting Yukawa interactions, this ratio is too large at the weak scale. For a top quark 
mass M t ~ 175 GeV, a good fit is obtained for small tan /3 ~ 1 [Dimopoulos et al (1992), 
Barger et al (1993)]. In this case, only the top quark Yukawa coupling is important. 
While for large tan/3 ~ 50 we recover the results of SO (10) Yukawa unification. For a 
recent analysis, see [Barr and Dorsner (2003)]. 
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Figure 11. Constant x 2 contours as a function of fi, Mx/2 for mie — 3 TeV. Note the 
much larger range of parameters with \ 2 < 1 for this larger value of vci\§. The green 
shaded region is consistent with the recent WMAP data for dark matter abundance of 
the neutralino LSP. The light shaded region in the lower left hand corner (separated 
by the solid line) is excluded by chargino mass limits, while the light shaded region in 
the upper left (right side) is excluded by a cosmological dark matter abundance which 
is too large (Higg mass which is too light). To the left of the vertical contour for a 
light Higgs with mass at the experimental lower limit, the light Higgs mass increases 
up to a maximum value of about 121 GeV at the lower left-hand acceptable boundary. 



5.2. Fermion mass hierarchy & Family symmetry 

In both the standard model and the MSSM, the observed pattern of fermion masses and 
mixing angles has its origin in the Higgs-quark and Higgs-lepton Yukawa couplings. In 
the standard model these complex 3x3 matrices are arbitrary parameters which are 
under constrained by the 13 experimental observables (9 charged fermion masses and 4 
quark mixing angles). [We consider neutrino masses and mixing angles in the following 
section.] In the MSSM more of the Yukawa parameters are in principle observable, 
since left and right-handed fermion mixing angles affect squark and slepton masses and 
mixing. [We consider this further in Section EH ] What can Bsay about fermion masses? 
Balone constrains the Yukawa sector of the theory simply by requiring that all terms in 
the superpotential are holomorphic. Combined with flavor symmetries, the structure of 
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fermion masses can be severely constrained. On the other hand, the only information we 
have about these flavor symmetries is the fermion masses and mixing angles themselves, 
as well as the multitude of constraints on flavor violating interactions. There are, 
perhaps, many different theories with different family symmetries that fit the precision 
low energy data (including fermion masses and mixing angles). The goal is to find a set 
of predictive theories, i.e with fewer arbitrary parameters than data, that fit this data. 
The more predictive the theory, the more testable it will be.^f Within the context of 
the MSSM, theories have been constructed with U(l) family symmetries [Binetruy et al 
(1996), Elwood et al (1997,1998), Irges et al (1998), Faraggi and Pati (1998), Kakizaki 
and Yamaguchi (2002), Dreiner et al (2003)], with discrete family symmetries [Frampton 
and Kephart (1995a,b), Hall and Murayama (1995), Carone et al (1996), Carone and 
Lebed (1999), Frampton and Rasin (2000), Aranda et al (2000)] or non-abelian family 
symmetries [Hall and Randall (1990), Dine et al (1993), Nir and Seiberg (1993), Pouliot 
and Seiberg (1993), Leurer et al (1993,1994), Pomarol and Tommasini (1995), Hall and 
Murayama (1995), Dudas et al (1995,1996), Barbieri et al (1996), Arkani-Hamed et al 
(1995,1996), Barbieri et al (1997), Eyal (1998)]. However, the most predictive theories 
combine both grand unified and family symmetries [Kaplan and Schmaltz (1994), Babu 
and Mohapatra (1995), Lucas and Raby (1996), Frampton and Kong (1996), Blazek 
et al (1997), Barbieri and Hall (1997), Barbieri et al (1997), Allanach et al (1997), 
Berezhiani (1998), Blazek et al (1999,2000), Derrmsek and Raby (2000), Shafi and 
Tavartkiladze (2000), Albright and Barr (2000,2001), Altarelli et al (2000), Babu et al 
(2000), Berezhiani and Rossi (2001), Kitano and Mimura (2001), Maekawa (2001), King 
and Ross (2003), Chen and Mahanthappa (2003), Raby (2003), Ross and Velasco-Sevilla 
(2003), Goh et al (2003), Aulakh et al (2003)]. The Yukawa couplings in a predictive 
theory are completely defined in terms of the states and symmetries of the theory. The 
ultimate goal of this program is to construct one (or more) of these predictive theories, 
providing good fits to the data, in terms of a more fundamental theory, such as M 
theory. Only then will higher order corrections to the theory be under full control. It 
is important to remark at this stage that any theory, derived from some fundamental 
theory, includes non-renormalizable higher dimension operators. The higher dimension 
operators are suppressed by the fundamental scale (for example, the string scale Ms) 
which is assumed to be greater than the GUT scale Mq. As we shall now see, these 
higher dimension operators are useful in explaining the hierarchy of fermion masses. 
The 3x3 up, down and charged lepton mass matrices are given by the mass terms: 

Crnass = u Y u u (H u ) + d Y d d (H d ) + e Y e e (H d ). (92) 

Empirical descriptions of the quark mass matrices have been discussed in all the papers 
referenced above. As an example, consider the following theory incorporating the 
hierarchy of masses and mixing angles in an SU(5) BGUT with U(l) family symmetry 

% Of course, if fior proton decay is observed then there will be much more low energy data available to 
test these theories. 
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Table 5. U(l) charge Q of Higgs and matter fields in the (1st, 2nd, 3rd) generation, 
field H u H d 10 = {Q, u, e} 5 = {d, L} V 
Q -2 1 (4,3,1) (4,2,2) (1,-1,0) 



by [Altarelli et al (2000)] 



( A 6 A 5 A 3 \ 

A 5 A 4 A 2 
\X 3 X 2 1 j 



Yd = Y e T 



( A 5 A 3 A 3 \ 

A 4 A 2 A 2 
\\> 1 1 ) 



A 4 . (93) 



Order 1 coefficients of the matrix elements are implicit. We then obtain the rough 
empirical relations. 

m c /m t ~ m u /m c ~ V 2 b ps A 4 (94) 
m s /m b ~ m d /m s ~ K u 2 ps A 2 . 

In addition, the Yukawa matrices for down quarks and charged leptons satisfy the SU(5) 
relations 

Afe = A T , X s = A M , Xd = X e . (95) 

This works for the third generation, as discussed in Section I5.1.3[ however it clearly 
doesn't work for the first and second generations, where it gives the unacceptable 
prediction 

20 ps m s /m d = m^/rrie ps 200. (96) 

Hence in this SU(5) model with U(l) family symmetry, an additional Higgs in the 
75 dimensional representation (with U(l) charge zero) also contributes to down quark 
and charged lepton Yukawa matrices [Altarelli et al (2000)]. The arbitrary, order one 
coefficients for each term in the Yukawa matrix are then fit to the quark masses and 
mixing angles and charged lepton masses. Note in this theory there are more arbitrary 
parameters, than fermion mass observables; hence there are no predictions for fermion 
masses and mixing angles. Nevertheless, predictions for proton decay are now obtained. 
Moreover, given a model for soft Bbreaking terms like the CMSSM, one can also predict 
rates for flavor violating processes. 

The structure for these Yukawa matrices are determined by the U(l) family 
symmetry, spontaneously broken by a scalar field <ft with U(l) charge —1. The symmetry 
breaking field is inserted into each element of the Yukawa matrix in order to obtain 
a U(l) invariant interaction. This results in effective higher dimensional operators 
suppressed by a scale M with A ~ ((f)) /M. This is the Froggatt-Nielsen mechanism 
[Froggatt and Nielsen (1979), Berezhiani (1983,1985), Dimopoulos (1983), Bagger et al 
(1984)]. For a review see [Raby (1995)]. Given the U(l) charge assignments in Table El 
[Altarelli et al (2000)] we obtain the Yukawa matrices in Eqn. IHB1 Similar mass matrices 
using analogous U(l) family symmetry arguments have also been considered. 
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Using a non-abelian family symmetry, such as SU(2) x U(l) or SU(3) or discrete 
subgroups of SU(2), models with fewer arbitrary parameters in the Yukawa sector have 
been constructed. An example of a very predictive SO(10) BGUT with SU(2) x U(l) ?l 
family symmetry is given by [Barbieri et al (1997(a,b),1999), Blazek et al (1999,2000)]. 
An analogous model can be obtained by replacing the SU(2) family symmetry with 
a discrete subgroup D(3) [Derrmsek and Raby (2000)]. The model incorporates the 
Froggatt-Nielsen mechanism with a hierarchy of symmetry breaking VEVs explaining 
the hierarchy of fermion masses. The effective fermion mass operators are given in Fig. 
[T2l In particular, the family symmetry breaking pattern 

SU 2 xUi — ► Ui — ► nothing (97) 

with small parameters e ~ e and e', respectively, gives the hierarchy of masses with 
the 3rd family 3> 2nd family 3> 1st family. It includes the Georgi - Jarlskog [Georgi 
and Jarlskog (1979)] solution to the unacceptable SU(5) relation (Eqn. I9ljj) with the 
improved relation 



m. 



(9£ 



(99) 



*n, m d ~ 3m e . 

This is obtained naturally using the VEV 

(45) = (B — L) M G . 

In addition, it gives the SO(10) relation for the third generation 

A t = \ b = \ T = \ Vr = A (100) 

and it uses symmetry arguments to explain why m u < m d even though m t 3> mj. 
Finally the SU2 family symmetry suppresses flavor violation such as p — > ej. When 
SO(10) x SU(2) x U(l) n is broken to the MSSM the effective Yukawa couplings (Eqn. 
I102j) are obtained. The superpotential for this simple model is given by 

W D 16 3 10 16 3 + 16 a 10 x a (101) 

ia cab 

+ Xa (M x X a + 45 %r 16 3 + 45 16 6 + A ab 16 6 ) 



where a , S ab = S ba , A ab = —A ba are the familon fields whose VEVs break the family 
symmetry, M x = M (1 + a X + (3 Y) with X, Y charges associated with U(l)x, y, the 
orthogonal U(l) subgroups of SO(10), and { x° \Xa } are the heavy Froggatt-Nielsen 
fields. After the heavy x states are integrated out of the theory we obtain the effective 
fermion mass operators given in Fig. These four Feynman diagrams lead to the 
following Yukawa matrices for quarks and charged leptons. 

/ e' p -e C \ 



Y u = 


-e' p e p -e 


A 




{ e£ e 1 


J 






( e' -ea i\ 




Y d = 


—e' e — e a 


A 




\*t e 1 J 





(102) 
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Xa 
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Figure 12. Effective fermion mass operators. The fields <?(> a , S"* 6 = S' ba ,v4 ab = -A ba 



spontaneously break the SU(2) x U(l) n family symmetry with e 
and e' oc (A 12 ). 



<5 2 



with 



/ 



\-3e£ -3 e 



3ee \ 



A 



> 1 )/(0 2 ); e oc (S 22 )/M; 



(103) 



e ex (0 2 )/M; 
1 + a 



1 — 3a' 



e' ~ (A i2 )/M; 
p ~ « a. 
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Table 6. Fit to fermion masses and mixing angles for SO(10) GUT with SU(2) x 
U(l) n family symmetry [Blazek et al (1999)]. 



Observable 


Data(er)masses in GeV 


Theory 


IvJ-Z 


91 1 87 (() 091 


91 1 7 


Mm 


80 388 (() 080") 


80 40 




1 1 664 (0 001 21 


1.166 


a -1 

"EM 


1 37 04 (0 141 


137.0 




fl 1 1 QD I'D 0031 


n 1 1 74 


n • 10 3 

l J new ±u 


-i 20 (^ 31 


+0.320 


M t 


173.8 (5.0) 


175.0 


m b (M b ) 


4.260 (0.11) 


4.328 


M b - M c 


3.400 (0.2) 


3.421 


m s 


0.180 (0.050) 


0.148 


m d /m s 


0.050 (0.015) 


0.0589 


Q- 2 


0.00203 (0.00020) 


0.00201 




1.777 (0.0018) 


1.776 




0.10566 (0.00011) 


.1057 


M e • 10 3 


0.5110 (0.00051) 


0.5110 


v us 


0.2205 (0.0026) 


0.2205 


v cb 


0.03920 (0.0030) 


0.0403 


v ub /v cb 


0.0800 (0.02) 


0.0691 


B K 


0.860 (0.08) 


0.8703 


B(b^ sj) -10 4 


3.000 (0.47) 


2.995 


TOTAL x 2 


3.39 





The model has only 9 arbitrary Yukawa parameters (6 real parameters 
{|A|, |e|, |e|, \p\, \a\, \e'\} and three phases {$ e = & P , *&a}) to fit the 13 fermion 
masses and mixing angles (we have taken £ = 0). The fit to the low energy data is 
given in Table More details of this fit are found in [Blazek et al (1999)] and the 
predictions for proton decay are found in [Derrmsek et al (2001)]. Note, the model fits 
most of the precision electroweak data quite well. In [Blazek et al (1999)] there is also 
a prediction for sin 2/3 = 0.54 which should be compared to the present experimental 
value 0.727 (0.036). The prediction for sin 2(3 is off by 5 o. In addition the present 
experimental value for V u b/V c b is 0.086 (0.008), hence this fit (Table is somewhat 
worse than before. Both of these quantities are predictions due solely to the zeros in 
the 11, 13 and 31 elements of the Yukawa matrices [Hall and Rasin (1993), Roberts 
et al (2001), Kim et al (2004)]. These poor fits are remedied with the addition of a 
non-vanishing 13 / 31 element, i.e. £ ^ 0. In this case a good fit is obtained with one 
additional real parameter [Kim et al (2004)]. 

+ Note, some of the data, used in this fit, have significantly improved in recent years. 
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Table 7. Input parameters at M G for fit in Table El [Blazek et al (1999)]. 



(l/a G , Mq, eg) = 
(A, e, cr, e, p, e', f ) = 
$ e = * e , $ P ) = 
(mie, M 1/2 , A , n(M z )) = 
((m Hd /m 16 ) 2 , (m Hu /m 16 ) 2 , tan/3) 



(24.52, 3.05 • 10 16 GeV, -4.08%), 
(0.79, 0.045, 0.84, 0.011, 0.043, 0.0031, 0.00), 
(0.73, -1.21, 3.72)rad, 

(1000, 300, -1437, 110) GeV, 
(2.22, 1.65, 53.7) 



5.3. Neutrino masses 

The combined data from all neutrino experiments can be fit by the hypothesis of neutrino 
oscillations with the neutrino masses and mixing angles given by 



Am 



atm 



m, 



m 



sin 20. 



m 



atm 



m-. 



3 x 10~ s eV z 



7 x 1(T 5 eV z 



(104) 



0.8 < sm26 so i < 1 

For recent theoretical analyses of the data, see [Barger et al (2003), Maltoni et al (2003), 
Gonzales-Garcia and Pena-Garay (2003)]. This so-called bi-large neutrino mixing is well 
described by the PMNS mixing matrix 

Csol Ssol \ / v\ \ 



V'2 
\ ^ ) 



s so i/V2 c sol /V2 l/y/2 
\ -s sol /V2 c so i/V2 -1/V2 J 
which takes mass eigenstates into flavor eigenstates. The 1-3 mixing angle satisfies 
sin# 13 < 0.2 at 3 a [Maltoni et al (2003)]. 

Using the See-Saw mechanism [Yanagida (1979), Glashow (1979), Gell-Mann et 
al (1979), Mohapatra and Senjanovic (1980)], neutrino masses are given in terms of 
two completely independent 3x3 mass matrices, i.e. the Dirac mass matrix m v and 
a Majorana mass matrix via the formula A4 l 



m,. 



M 



N 



m,, 



The smallness of 

neutrino masses is explained by the large Majorana mass scale, of order 10 14 — 10 15 GeV; 
very close to the GUT scale. In addition, the large mixing angles needed to diagonalize 
}A V can be directly related to large mixing in m„, in or in some combination of 
both. Lastly, the Dirac neutrino mass matrix m u is constrained by charged fermion 
masses in SO(10), but not in SU(5) where it is completely independent. 

The major challenge for theories of neutrino masses is to obtain two large mixing 
angles; as compared to charged fermions where we only have small mixing angles in 
Vckm- There are several interesting suggestions for obtaining large mixing angles in the 
literature. [For recent reviews of models of neutrino masses, see [Altarelli and Feruglio 
(2003), Altarelli et al (2003), King (2003)].] 

• Degenerate neutrinos and RG running [Mohapatra et al (2003), Casas et al (2003)], 
It was shown that starting with three degenerate Majorana neutrinos and small 
mixing angles at a GUT scale, that RG running can lead to bi-large neutrino 
mixing at low energies. 
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• Minimal renormalizable SO(10) [Goh et al (2003), Bajc et al (2003)], 

SO(10) with Higgs in the 10 and 126 representations can give predictable theories 
of fermion masses with naturally large neutrino mixing angles. 

• Dominant Majorana neutrinos [King (1998,2000)], 

It was shown that large neutrino mixing can be obtained via coupling to a single 
dominant right handed neutrino. 

• Minimal Majorana sector [Frampton et al (2002), Raidal and Strumia (2003), Raby 
(2003)], 

It was shown that a simple model with two right handed neutrinos can accommodate 
bi-large neutrino mixing with only one CP violating phase. In such a theory, CP 
violating neutrino oscillations measured in low energy accelerator experiments are 
correlated with the matter - anti-matter asymmetry obtained via leptogenesis. 

• Lopsided charged lepton and down quark matrices [Lola and Ross (1999), Nomura 
and Yanagida (1999), Albright and Barr (2000 (a,b), 2001), Altarelli et al (2000), 
Barr and Dorsner (2003)]. 

In SU(5) (or even in some SO(10) models) the down quark mass matrix is related 
to the transpose of the charged lepton mass matrix. A large left-handed fi — r 
mixing angle is thus directly related to a large right-handed s — b mixing angle. 
Whereas right-handed quark mixing angles are not relevant for CKM mixing, the 
large left-handed charged lepton mixing angle can give large — v T mixing. 

Let us now consider the last two mechanisms in more detail. 



5.3.1. SU(5) x U(l) flavor symmetry One popular possibility has the large — v T 
mixing in the Dirac charged Yukawa matrix with 

/ A 5 A 4 A 2 \ 



Y P 



A 3 
A 3 



A 2 
A 2 



(105) 



The neutrino Dirac Yukawa matrix and the Majorana matrix are given by 
/A 3 A A 2 \ 



V 



M, 



N 



A 
A 

( A 2 
1 








1 A\ 


1 







(106) 



(107) 



where we use the results of [Altarelli et al (2000)]. The light neutrino mass matrix is 
given by the standard See-Saw formula. We obtain: 

/ A 4 A 2 A 2 \ 

U e v 2 JM. (108) 



M„ = U. 



tr 



A 2 
A 2 
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where v u is the VEV of the Higgs doublet giving mass to the up quarks, M is the 
heavy Majorana mass scale and all entries in each matrix are specified up to order one 
coefficients. U e is the mixing matrix taking left-handed charged leptons into the mass 
eigenstate basis. It is given by 



/ 












m. 






mi 



where 



( A 4 

A 2 

v^ 2 



\ 

J 

A 2 
1 
1 



U? (Y e r e t) U* e (H d ) 



(109) 



A 2 \ 
1 

1 J 



(110) 



up to order one coefficients. It is clear that without order one coefficients, the neutrino 
mixing matrix is the identity matrix. Hence the arbitrary order one coefficients are 
absolutely necessary to obtain bi-large neutrino mixing.* Hierarchical neutrino masses 
with m 3 ^> m 2 m\ and bi-large neutrino mixing can naturally be obtained [Altarelli 
and Feruglio (2003)]. 

Of course, different U(l) charge assignments for all the fields can lead to other 
experimentally acceptable solutions to the solar neutrino problem. For a review and 
further references, see [Altarelli and Feruglio (2003)]. 



5.3.2. SO 10 x[SU 2 xU?] FS model Within the context of the SO 10 x[SU 2 xU?] FS model, 
bi-large neutrino mixing is naturally obtained using the mechanism of the minimal two 
Majorana neutrino sector. The Dirac neutrino mass is fixed once charged fermion masses 
and mixing angles are fit. It is given by the formula: 



Y„ 







e' uj 



V -3 e f a 



-e' UJ 
3 e uj 
— 3 t o~ 



e uj 



A 



with u = 2 <t/(2 cr — 1) and the Dirac neutrino mass matrix given by m u = Y u -^ sin (5. 
Of course, all the freedom is in the Majorana neutrino sector. The FGY ansatz 
[Frampton et al (2002)] is obtained with the following Majorana neutrino sector [Raby 
(2003)]: 



16 



neutrino 



r (N, r 16 Q + N 2 r 16 a + N 3 9 16 3 ) 
M v ' 

+ l - (S, Nl + S 2 Nj) 

where {A^, % = 1,2,3} are SO(10) and SU(2) - flavor singlets. In this version of the 
theory, the symmetric two index tensor flavon field S ab is replaced by an SU(2) doublet 
4> a such that S ab = 4> a 4> b / ' M . Note, since the singlet N$ has no large Majorana mass, it 
gets a large Dirac mass by mixing directly with z/ 3 at the GUT scale. Thus u 3 is removed 

* I thank G. Altarelli, private communication, for emphasizing this point. 
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from the See-Saw mechanism and we effectively have only two right-handed neutrinos 
taking part. 

Integrating out the heavy neutrinos we obtain the light neutrino mass matrix given 

by 

M v = U l : \D tr M N X D\XJ e (111) 
where U e is the unitary matrix diagonalizing the charge lepton mass matrix and 

/ a \ 

, Mn = 



D 



tr 



a' b 
b' 



(Si) 
(S 2 ) 



with 



/ , / , r , , , 1 v \ M v sin 3 

b= € UJ X (M 2 / (0 )) jf 

, , , i . . M v sin 3 
b >= -3e^A(M 2 /(0 1 )) 7 f. 



e' u A (M!/(0 2 )) 



2 MvsinS 



Vie y/2 



a' = (-e' r 1 + 3£>A (M 1 /(0 2 )) 



M v sin /3 

fl6 V 7 ^ 



We obtain 

b ~ 6' 

naturally since e' ~ e £. In addition we can accommodate 
a ~ a' 



(112) 



(113) 



(114) 



115) 



with minor fine-tuning 0(1/10) since e' C,^ 1 ~ e. Note, this is the Frampton-Glashow- 
Yanagida ansatz [Frampton et al (2002)] with a bi-large neutrino mixing matrix obtained 
naturally in a BGUT. 



5.4- Flavor Violation 

Quarks and leptons come in different flavors: up, down, charm, strange, top, bottom; 
electron, muon, tau. We observe processes where bottom quarks can decay into charm 
quarks or up quarks. Hence quark flavors (for quarks with the same electric charge) 
are interchangeable. In the standard model, this is parametrized by the CKM mixing 
matrix. We now have direct evidence from neutrino oscillation experiments showing 
that tau, muon and electron numbers are not separately conserved. Yet, we have never 
observed muons changing into electrons. In supersymmetric theories there are many 
more possible ways in which both lepton and quark flavors can change. This is because 
scalar quarks and leptons carry the flavor quantum numbers of their fipartners. Thus 
flavor violation in the scalar sector can lead to flavor violation in the observed fermionic 
sector of the theory. This gives rise to the SUSY flavor problem. We consider two 
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Figure 13. The one loop contribution to the process \i — > proportional to an 
off-diagonal scalar muon - electron mass term in the charged lepton mass eigenstate 
basis. 



Table 8. Some constraints from the non-observation of flavor violation on squark, 
slepton and gaugino masses [Gabbiani et al (1996)]. For the electron electric dipole 
moment we use the relation d e N ~ 2(100/m ; -(GeV)) 2 sm$ J 4 i B x 10 _23 e cm. 

Observable Experimental bound (i) Experimental bound (2) 

B(/ 1 ^ e7 )<1.2xl0- 11 \{5[ 2 )ll\ < 2.1 x 10- 3 (m r (GeV)/f00) 2 \{5[ 2 )ll \ < 0.8 (m r (TcV)/2) 2 

Am K < 3.5 x f(T 12 MeV J \Re(5f 2 ) 2 LL \ < 1 .9 x fO- 2 (m,j(GeV)/500) \j \Re{5f 2 ) 2 LL \ < 7.6 x 10~ 2 (m§(TeV)/2) 

e K < 2.28 x lO- 3 \J\Im(Sf 2 ) 2 LL \ < 1 .5 x fO- 3 (m,(GcV)/500) \J \Re{5f 2 ) 2 LL \ < 6.0 x 10~ 3 (m ? (TeV)/2) 

d e N < 4.3 x f0" 27 e cm sin$ A ,B < 4 x f0~ 4 x (m f (GeV)/100) 2 sin$ A , B < 0.16 x (to,- (TeV)/2) 2 



examples here: fi — > or B s — > We show why SUSY GUTs and/or neutrino 

masses can cause enhanced flavor violation beyond that of the standard model. In this 
section we consider several ways to solve the SUSY flavor problem. 

But first let us illustrate the problem with a few of the dominant examples [Gabbiani 
et al (1996)]. In the first column of Table |H1 we present four flavor violating observables 
with their experimental bounds. In the second and third columns we present the bounds 
on flavor violating scalar mass corrections 

8(j = Am^/m 2 (116) 

where i, j = 1, 2, 3 are family indices, Am?- is an off-diagonal scalar mass insertion for 
/ = {quark, lepton} flavor (treated to lowest non-trivial order in perturbation theory) 
in the flavor diagonal basis for quarks and leptons. fa is the average squark or slepton 
mass squared. The subscripts LL refer to left-handed squark or slepton mass insertions. 
There are separate limits on RR and LR mass insertions [Gabbiani et al (1996)] which 
are not presented here. The only difference between the second and third columns is 
the fiducial value of m 2 . Clearly as the mean squark or slepton mass increases, the fine 
tuning necessary to avoid significant flavor violation is dramatically reduced. As seen 
in Fig. the amplitude for \i — > is proportional to 5f 2 and is suppressed by 1/m 2 , 
since it is an effective dimension 5 operator. 
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5.4-1- The Origin of flavor violation in fitheories There are three possible ways to 
avoid large flavor violation. 

(i) Having squarks and sleptons, with the same standard model gauge charges, be 
degenerate and, in addition, the cubic scalar interactions proportional to the 
Yukawa matrices. 

(ii) Alignment of squark and slepton masses with quark and lepton masses. 

The fermion and scalar mass matrices are "aligned" when, in the basis where 
fermion masses are diagonal, the scalar mass matrices and cubic scalar interactions 
are approximately diagonal as well. 

(iii) Heavy first and second generation squarks and sleptons. 

The CMSSM (or mSUGRA) is an example of the first case. It has a universal scalar 
mass itlq and tri-linear scalar interactions proportional to Yukawa matrices. These 
initial conditions correspond to a symmetry limit [Hall et al (1986)] - dubbed minimal 
flavor violation [Ciuchini et al (1998)]- where the only flavor violation occurs in the 
CKM matrix at the messenger scale for Bbreaking, or in this case, the Planck scale. 
Gauge-mediated Bbreaking, where squarks and sleptons obtain soft Bbreaking masses 
via standard model gauge interactions, is another example of the first case (for a 
review, see [Giudice and Rattazzi (1999)]. In this case the messenger mass is arbitrary. 
Finally, within the context of perturbative heterotic string theory, dilaton Bbreaking 
gives universal scalar masses at the string scale. For moduli Bbreaking, on the other 
hand, scalar masses depend on modular weights, whose values are very model dependent. 

Abelian flavor symmetries can be used to align quark (lepton) and the corresponding 
squark (slepton) mass matrices, but they still require one of the above mechanisms for 
obtaining degenerate scalar masses at zeroth order in symmetry breaking. Non-abelian 
symmetries, on the other hand, can both align fermion and scalar mass matrices and 
guarantee the degeneracy of the scalar masses at zeroth order. 

Finally, since the most stringent limits from flavor violating processes come from 
the lightest two families, if the associated squarks and sleptons are heavy these processes 
are suppressed [Dimopoulos and Giudice (1995)]. A natural mechanism for obtaining 
this inverted scalar mass hierarchy with the first and second generation scalars heavier 
than the third was discussed by [Bagger et al (1999,2000)]. 

If the messenger scale for Bbreaking is above the GUT scale or even above the See- 
Saw scale for neutrino masses then squark and slepton masses can receive significant 
flavor violating radiative corrections due to this beyond the standard model physics 
[Hall et al (1986), Georgi (1986), Borzumati and Masiero (1986), Leontaris et al (1986), 
Barbieri et al (1995a,b), Hisano et al (1995,1996)]. In addition, if the fermion mass 
hierarchy is due to flavor symmetry breaking using a Froggatt-Nielsen mechanism, 
then upon integrating out the heavy Froggatt-Nielsen sector new flavor violating soft 
Bbreaking terms may be induced [Dimopoulos and Pomarol (1995), Pomarol and 
Dimopoulos (1995)]. Hence the low energy MSSM is sensitive to physics at short 
distances. This is both a problem requiring natural solutions and a virtue leading 
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to new experimentally testable manifestations of Btheories. 

For example in BGUTs, color triplet Higgs fields couple quarks to leptons. As 
a consequence flavor mixing in the quark sector can, via loops, cause flavor mixing, 
proportional to up quark Yukawa couplings, in the lepton sector. This comes via off- 
diagonal scalar lepton masses in the basis where charged lepton masses are diagonal. 
While the one loop contribution of charm quarks give a branching ratio Br(fi — > ej) ~ 
1CT 15 [Hall et al (1986)], the top quark contribution leads to very observable rates 
[Barbieri et al (1995a)] near the experimental bounds. Moreover, new experiments 
will soon test these results. In addition, experimental evidence for neutrino oscillations 
makes it clear that the lepton sector has its own intrinsic flavor violation. In the standard 
model, these effects are suppressed by extremely small (< eV) neutrino masses. In 
Bhowever, flavor violation in the (s) neutrino sector leads, again via loops, to eminently 
observable mixing in the charged (s)lepton sector [Hisano et al (1995,1996)]. There 
are a large number of papers in the literature which try to use low energy neutrino 
oscillation data in an attempt to predict rates for lepton flavor violation. However 
a bottom-up approach is fraught with the problem that low energy oscillation data 
cannot completely constrain the neutrino sector [Casas and Ibarra (2001), Lavignac 
et al (2001,2002)]. It has been shown that neutrino oscillation data and U — > lj 7 
measurements can nevertheless provide complementary information on the See-Saw 
parameter space [Davidson and Ibarra (2001), Ellis et al (2002)]. In a recent analysis, it 
was shown that lepton flavor violation can constrain typical BSO(10) theories [Masiero 
et al (2003)]. Finally, it is important to note, that the same physics can lead to enhanced 
contributions to flavor conserving amplitudes such as the anomalous magnetic moment 
of the neutrino (a M ) [Chattopadhyay and Nath (1996), Moroi (1996)] and the electric 
dipole moments of the electron (df) and neutron (d™) [Dimopoulos and Hall (1995), 
Hisano and Tobe (2001), Demir et al (2003)]. Moreover, the rates for these flavor 
violating processes increase with tan f3. 

Let us now consider flavor violating hadronic interactions at large tan/5. We 
focus on a few important examples, in particular, the processes B — > X s ^f, B — > 
X s l + l~ forward-backward asymmetry and B s — > For a more comprehensive 

study, see [Hall et al (1994), Hempfling (1994), Carena et al (1994), Blazek et al 
(1995), Chankowski and Pokorski (1997), Misiak et al (1998), Huang and Yan (1998), 
Huang et al (1999), Hamzaoui et al (1999), Babu and Kolda (2000), Chankowski 
and Slawianowska (2001), Carena et al (2001), Bobeth et al (2001), Huang et al 
(2001), Dedes et al (2001), Isidori and Retico (2001), Buras et al (2002,2003), 
Dedes and Pilaftsis (2003)]. The standard model contribution to B — > X s 7 has 
significant uncertainties, but the calculated branching ratio is consistent with the 
latest experimental data. Supersymmetry contributions are typically divided into two 
categories, i.e. the contributions contained in a two Higgs doublet model and then 
the rest of the Bspectrum. The charged Higgs contribution has the same sign as the 
standard model contribution and thus increases the predicted value for the branching 
ratio. This spoils the agreement with the data and thus a lower limit on the charged 
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Figure 14. The one loop chargino contribution to the process B — > X s j proportional 
to At and tan/3. 



Higgs mass is obtained. In the minimal flavor Bscenario, the additional Bcontribution is 
dominated by the chargino loop (Fig. The sign of this term depends on the sign 
of \i. For n > (this defines my conventions) the chargino contribution is the opposite 
sign of the standard model contribution to the coefficient Cj of the magnetic moment 
operator Oj ~ sl E^j, or . Moreover, this contribution is proportional to tan/?. 
For small or moderate values of tan f3 the ficorrection to CV is small and for fi > it 
tends to cancel the charged Higgs and standard model contributions. This is in the right 
direction, giving good agreement with the data. For [i < the agreement with the data 
gets worse and can only work for large Higgs and squark masses, so that the overall 
Bcontribution is small. On the other hand, for fi > and large tan/3 ~ 50 there is 
another possible solution with the total Bcontribution to CV equal to twice the standard 
model contribution but with opposite sign. In this case C\ otal = Cf +Cj ~ —C® and 
good fits to the data are obtained [Blazek and Raby (1999)]. Although the sign of CV is 
not observable in B s — > X s 7, it can be observed by measuring the forward - backward 
asymmetry in the process B — > X s l + l~ [Huang and Yan (1998), Huang et al (1999), 
Lunghi et al (2000), Ali et al (2002), Bobeth et al (2003)] where forward (backward) 
refers to the positive lepton direction with respect to the B flight direction in the rest 
frame of the di-lepton system. f 

5.4-2. Flavor violating Higgs couplings at large tan/3 The MSSM has two Higgs 
doublets which at tree level satisfy the Glashow - Weinberg condition for natural flavor 
conservation. Up quarks get mass from H u and down quarks and charged leptons get 
mass from Hj. Thus when the fermion mass matrices are diagonalized (and neglecting 
small neutrino masses) the Higgs couplings to quarks and leptons are also diagonal. 
However, this is no longer true once Bis broken and radiative corrections are considered. 
In particular, for large values of tan/3 the coupling of H u to down quarks, via one 
loop corrections, results in significant flavor violating vertices for neutral and charged 

f I thank K. Tobe for pointing out this possibility to me. 
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Higgs. These one loop corrections to the Higgs couplings contribute to an effective 
Lagrangian (Eqn. I117j) [Blazek et al (1995), Chankowski and Pokorski (1997)]. The 
chargino contribution (Fig. EJ) is proportional to the square of the up quark Yukawa 
matrix, which is not diagonal in the diagonal down quark mass basis. As a result, at 
one loop order, the down quark mass matrix is no longer diagonal (Eqn. I118j) . This 
leads to tan f3 enhanced corrections to down quark masses and to CKM matrix elements 
[Blazek et al (1995)]. Upon re-diagonalizing the down quark mass matrix we obtain 
the effective flavor violating Higgs - down quark Yukawa couplings given in Eqns. 11191 
and 11201 [Chankowski and Pokorski (1997), Babu and Kolda (2000), Chankowski and 
Slawianowska (2001), Bobeth et al (2001), Huang et al (2001), Dedes et al (2001), 
Isidori and Retico (2001), Buras et al (2002,2003), Dedes and Pilaftsis (2003)]. 

du Xt 9 d Rl H° d * (117) 



rddH 



- d Li A \j d Rj H d 



0* 



du 8X1 d Rj H° + h.c. 



m 



Diagonal 



V 



+ AA d + <5A d tan/3 



v cos (3 



'1181 



where 



- T2 3 '- 



Ft Pr + f>;; p. 



F " p R + F " p i. 



li* 



F A P R + F A * Pr 



d' 3 H 



d) A, 



(119) 



1 + tan 2 f3) cos/? cos(a — (3), 
1 + tan 2 (3) cos/3 sin(a — (3), 



i.i 



<JAj'(l + tan 2 /3) cos (3. 



(120) 



This leads to tan (3 enhanced flavor violating couplings for the neutral Higgs bosons. 
For example, the branching ratio B(B S — ► /U + /i~) is proportional to tan (3 A and inversely 
proportional to the fourth power of the CP odd Higgs mass uia [Babu and Kolda (2000), 
Chankowski and Slawianowska (2001), Bobeth et al (2001), Huang et al (2001), Dedes 
et al (2001), Isidori and Retico (2001), Buras et al (2002,2003), Dedes and Pilaftsis 
(2003)], since the contributions of the two CP even Higgs bosons approximately cancel 
[Babu and Kolda (2000)]. The present DO and CDF bounds constrain rriA > 250 
GeV for tan/3 ~ 50, although this result is somewhat model dependent [Dermfsek 
et al (2003)]. Note the DO and CDF bounds from the Tevatron Run 2 now give 
B(B S -> fj+fj,-) < 1.2 x 10~ 6 (CDF) [Lin (2003)] at 95% CL and < 1.6 x 10~ 6 (DO) 
[Kehoe (2003)]. An order of magnitude improvement will test large tan/3 fifor up to 
~ 500 GeV [Dermfsek et al (2003)] (see Fig. [TJ. 
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m ia =3TeV, m A =500 GeV 
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Figure 15. The dashed (blue) lines are contours for constant branching ratio 
B(B S — > /i + /i _ ) as a function of /i, Mi/ 2 for mi 6 = 3 TeV. The green shaded region is 
consistent with the recent WMAP data for dark matter abundance of the neutralino 
LSP. The light shaded region in the lower left hand corner (separated by the solid line) 
is excluded by chargino mass limits, while the light shaded region in the upper left 
(right side) is excluded by a cosmological dark matter abundance which is too large 
(Higg mass which is too light). 

6. fiDark Matter 

The two most popular dark matter candidates are axions or the LSP of SUSY. Both are 
well motivated cold dark matter candidates. There have been several recent studies of 
SUSY dark matter in light of the recent WMAP data. For these analyses and references 
to earlier works, see [Roszkowski et al (2003), Ellis et al (2003c, d,e,f), Chattopadhyay 
et al (2003)]. These calculations have been performed with different assumptions about 
soft Bbreaking parameters, assuming the CMSSM boundary conditions at Mq [Ellis et al 
(2003bc,d,f), Chattopadhyay et al (2003)] or arbitrary low energy scalar masses [Ellis et 
al (2003d)]. Soft Bbreaking outside the realm of the CMSSM has also been considered. 
For example, soft breaking with non-universal Higgs masses have been analyzed recently 
by [Ellis et al (2003a,b), Roszkowski et al (2003)]. In the latter case, the soft Bbreaking 
parameters consistent with SO(10) Yukawa unification were studied. 
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In the limit of large squark and slepton masses, it is important to have efficient 
mechanisms for dark matter annihilation. Most recent studies have focused on the 
neutralino LSP in the limit of large tan/3 and/or the focus point limit. In both cases 
there are new mechanisms for efficient neutralino annihilation. For large tan/5 > 40 
neutralino annihilation via direct s-channel neutral Higgs boson exchange dominates 
[Roszkowski et al (2001), Ellis et al (2001a,b)]. In this limit the CP even and odd Higgs 
bosons have large widths, due to their larger coupling to bottom quarks and r leptons. 
In the focus point limit, on the other hand, the neutralino LSP is a mixed Higgsino - 
gaugino state. Thus it has more annihilation channels than the pure bino LSP case, 
valid for values of the universal scalar mass mo < TeV [Feng et al (2000c, 2001)]. 

Direct detection [Ellis et al (2003b,e), Roszkowski et al (2003), Munoz (2003)] 
and/or indirect detection [Baer and Farrill (2003), de Boer et al (2003)] of neutralino 
dark matter has also been considered. In fact, [de Boer et al (2003)] suggests that some 
indirect evidence for Bdark matter already exists. 

7. Open questions 

It is beyond the scope of this review to comment on many other interesting topics affected 
by supersymmetric theories. Several effective mechanisms for generating the matter- 
anti-matter asymmetry of the universe have been suggested, including the Affleck-Dine 
mechanism [Affleck and Dine (1985)], which is purely a supersymmetric solution, or 
leptogenesis [Fukugita and Yanagida (1986)], which is not necessarily supersymmetric. 
There have also been many interesting studies of inflation in a Bcontext. Finally, we 
have only made passing reference to superstring theories and Bbreaking mechanisms or 
fermion masses there. 

Simple " naturalness" arguments would lead one to believe that SUSY should have 
already been observed. On the other hand, "focus point" or "minimal SO(10) SUSY" 
regions of soft Bbreaking parameter space extend to significantly heavier squark and 
slepton masses without giving up on "naturalness. " In both the "focus point" and 
"mSOioSM" regions of parameter space we expect a light Higgs with mass of order 
114 — 120 GeV. Both ameliorate the Bflavor problem with heavy squark and slepton 
masses. They are nevertheless both surprisingly consistent with cosmological dark 
matter abundances. In addition, we have shown that the "mSOi SM" satisfies Yukawa 
coupling unification with an inverted scalar mass hierarchy. Thus one finds first and 
second generation squarks and sleptons with mass of order several TeV, while gauginos 
and third generation squarks and sleptons are much ligher. In addition, it requires large 
values of tan /3 ~ 50 resulting in enhanced flavor violation. 

SUSY GUTs are the most natural extensions of the standard model, and thus 
they are the new "standard model" of particle physics. The mSUGRA (or CMSSM) 
boundary conditions at the GUT scale provide excellent fits to precision low energy 
electroweak data. BGUTs, besides predicting gauge coupling unification, also provide a 
framework for resolving the gauge hierarchy problem and understanding fermion masses 
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and mixing angles, including neutrinos. It also gives a natural dark matter candidate, 
and a framework for leptogenesis and inflation. 

BUT there are two major challenges with any supersymmetric theory. We do not 
know how Bis spontaneously broken or the origin of the // term. We are thus unable 
to predict the Bparticle spectrum, which makes Bsearches very difficult. Nevertheless, 
"naturalness" arguments always lead to some light Bsector, observable at the LHC, 
a light Higgs, with mass less than 0(135 GeV), or observable flavor violating rates 
beyond that of the standard model. Assuming Bparticles are observed at the LHC, 
then the fun has just begun. It will take many years to prove that it is really 
supersymmtry. Assuming Bis established, a Bdesert from Mz to Mq (or Mjv) becomes 
highly likely. Thus precision measurements at the LHC or a Linear Collider will probe 
the boundary conditions at the very largest and fundamental scales of nature. With the 
additional observation of proton decay and/or precise GUT relations for sparticle masses, 
SUSY GUTs can be confirmed. Hence with experiments at TeV scale accelerators or 
in underground detectors for proton decay, neutrino oscillations or dark matter, the 
fundamental superstring physics can be probed. Perhaps then we may finally understand 
who ordered three families. It is thus no wonder why the elementary particle physics 
community is desperately seeking SUSY. 
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